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We present a new chiral expansion scheme for the nucleon-nucleon scattering am- 
plitude which preserves unitarity exactly. Our effective field theory builds on the 

pg . power counting rules for 2-nucleon reducible diagrams proposed in [1]. We evalu- 

ate the leading order terms of the isospin one scattering amplitude and elaborate 
in detail on the 1 5o phase shift. Our chiral description of the 1 S'o-phase shift does 
compete in quality with modern phenomenological nucleon-nucleon potentials. We 
describe elastic and inelastic scattering quantitatively up to laboratory energies of 

g ; £ lab ~ 600 MeV. 
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q ' 1 Introduction 

The nucleon-nucleon scattering problem has a long history in nuclear physics 
/\ [2]. Nowadays there exist phenomenological nucleon-nucleon potentials [3-6] 

which successfully describe scattering data to high precision. These potentials 
are commonly constructed applying meson-exchange phenomenology. A seri- 
ous drawback of the traditional approach is, however, the lack of a systematic 
scheme in the sense that there does not exist a well defined procedure to im- 
prove the potential. One has not yet identified a small expansion parameter 
which guides the construction of the potential nor has one resolved the renor- 
malization problem. Results depend explicitly on ad-hoc formfactors. This will 
be an interesting issue when the data base on nucleon-nucleon scattering is 
further improved in the near future, in particular with upcoming experiments 
at the proton cooler facility COSY in Jiilich. 

A first attempt to systematically solve the problem of nucleon-nucleon scat- 
tering and construct a bridge to QCD was made by Weinberg [7]. He sug- 
gested to derive a nucleon-nucleon potential in time-ordered chiral perturba- 
tion theory [7]. The heavy meson-exchange contributions are, in accordance 
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with chiral symmetry, systematically absorbed in local two-nucleon interaction 
vertices. This program was first carried out by Ordonez, Ray and van Kolck 
[8,9] who computed nucleon phase shifts from the solution of a cutoff regu- 
larized Schrodinger equation feeded properly with the chiral potential. More 
recent attempts within the potential approach can be found in [10-12]. The 
chiral potential scheme is, however, still plagued with serious shortcomings. 
The systematic renormalization of the chiral potential scheme remains an open 
problem. This leads again to a regularization scheme dependence of physical 
results. Though in certain cases the regularization scheme dependence is found 
to be weak the reason for which awaits a systematic explanation. 

In [1] the author suggested to apply a relativistic version of chiral perturbation 
theory (%PT) with power counting rules also for two nucleon reducible dia- 
grams. The proposed scheme is practicable since the 'difficult' non-local pion 
dynamics is found to be perturbative. Non perturbative effects are generated 
by properly renormalized local two-nucleon vertices. As a result the s-wave 
scattering amplitude is the sum of a pole term ( representing the deuteron 
bound state in the spin triplet channel or the pseudo-bound state in the spin 
singlet channel ) and a smooth remainder which is evaluated perturbatively 
in the pion dynamics. The question whether pions are perturbative in the nu- 
cleon scattering problem has been addressed also in [13,14]. Our scheme in 
[1] shows similarities with a framework proposed recently by Kaplan, Savage 
and Wise [15] and questioned by Cohen and Hansen [16]. Starting with the 
static one-pion exchange potential the authors [15] arrive at the power count- 
ing rules for reducible diagrams suggested in [1]. However, as will be shown 
in detail the implementation of the generalized chiral counting rules of [1] 
in [15] is in contradiction to the low energy structure of the nucleon-nucleon 
scattering amplitude. In particular unitarity is strongly violated. The s-wave 
scattering amplitude acquires a more complicated structure. Its representation 
in terms of a pole term and a smooth remainder accessed in terms of perturba- 
tive pions fails at unexpectedly small momenta. We will demonstrate that the 
failure of this scheme [1] in the spin singlet channel is caused by a second non- 
perturbative phenomenon present in the subthreshold scattering amplitude. 
The purpose of this work is to develop a consistent effective field theoretic 
application of the chiral power counting rules proposed in [1] in accordance 
with the low energy structure of the scattering amplitude. 

We choose the relativistic version of the chiral Lagrangian since it offers on the 
one hand a manifest covariant framework and on the other hand it is the natu- 
ral scheme to describe particle production processes like pp — > pp7i°,pnn + . A 
further important novel ingredient of our scheme is a reorganization of the lo- 
cal two-body interaction vertices. Rather than expanding the interaction terms 
around zero three-momenta we propose to expand around a finite value. The 
physical motivation for this reorganization follows from the empirical obser- 
vation that for example the 1 Sq partial wave shows a zero at intermediate 



energies Ei a b. ~ 280 MeV. A novel systematic approximation scheme for the 
solution of the Bethe-Salpeter equation consistent with the chiral counting 
rules and unitarity is presented. 

Our work is subdivided into four parts. In section two we discuss the chiral 
counting rules as they emerge from the relativistic chiral Lagrangian and intro- 
duce our expansion scheme for the nucleon-nucleon scattering amplitude. In 
section three we present analytic results of explicit evaluations at subleading 
orders. By evaluating the old scheme of [1] at three-loop level it is demon- 
strated that indeed it fails at momenta p > m n . Many technical details of this 
calculation can be found in a series of appendices. In section four we extend 
our scheme to describe the pion production process and confront our results 
with the empirical scattering phase shift. We obtain a quantitative description 
of elastic and inelastic s-wave nucleon-nucleon scattering in the spin singlet 
channel up to E iab ~ 600 MeV. 



2 Discussion of relativistic chiral power counting rules 



The chiral Lagrangian is constructed conveniently in terms of Weinberg's stere- 
ographic coordinates [21]. The buildings blocks are the pionic four vector D^ 
with 



D c _ «^7T C 



" 1 + A% 2 

and the relativistic nucleon field N. The coupling A is identified with the chiral 
limit value of the pion decay constant A -1 = 2 f n ~ 186 MeV. The chiral 
Lagrangian is now the infinite sum of all isoscalar terms formed with any 
number of the pion field D^, nucleon fields N and their respective covariant 
derivatives with: 



V»N = [d, + i Y^- 2 r ■ (vf x d^ N , 

V,D a u = dpDl + 2 \n a (A, • 4) - 2 A (V • D v ) D^ . (2) 

Further terms are to be added in order to describe the explicit chiral symme- 
try breaking pattern due to small current quark masses. Such terms do not 
add further complications to our discussion and therefore are not considered 
explicitly here. Weinberg's chiral power counting rules [7,21] induce a system- 
atic truncation of this infinite hierarchy. At a given order only a finite number 
of Feynman diagrams derived from a finite number of interaction terms need 



to be considered. Since we will generalize Weinberg's counting rules we briefly 
recall his scheme. 

For a given irreducible Feynman diagram one introduces the number of internal 
nucleon lines, In, the number of internal pion lines, /„., the number of loops, 
L and the number, Vi, of vertices of type % with di 'small' derivatives involved. 
Here one calls a derivative small if it acts on the pion field. The chiral counting 
rule for a given Feynman diagram results if the nucleon propagator is given 
the chiral power —1, the pion propagator —2, a small derivative the power 

1 and the loop momentum the power 4. The usefulness of the counting rules 
(3) follows by means of the topological relations L = 1^ + In — J2i Vi + 1 and 

2 In + En = J2Vi n i where the number of external nucleon lines En and rii, 
the number of nucleon fields in an interaction term of type % is introduced. 
The well know result 



v = 4L-I N -2I 7T + J2V t d t 

i 

= 2- 1 1 E N + 2L + Y,V(d i + 1 1 n l -2) (3) 

i 

shows that for di + |rij — 2 > 0, higher loop diagrams are more and more 
suppressed since they carry larger and larger chiral powers. We emphasize 
that the counting rule (3) is not applicable for multi-nucleon reducible dia- 
grams which one encounters for example in the iteration of the Bethe-Salpeter 
equation for the nucleon-nucleon scattering amplitude. 

There are two problems inherent with the relativistic approach. First, any 
covariant derivative acting on the nucleon field produces the large nucleon 
mass and therefore must be assigned the minimal chiral power zero. Thus an 
infinite tower of interaction terms needs to be evaluated at given finite chiral 
order [17,18]. Second, the straightforward evaluation of Feynman diagrams in- 
volving relativistic nucleon propagators generates positive powers of the large 
nucleon mass [17]. Therefore the chiral power counting rule (3) is spoiled. The 
relativistic scheme appeared impractical and the heavy mass formulation of 
XPT, which overcomes both problems by performing a non relativistic 1/m 
expansion at the level of the Lagrangian density, was developed [19] and ap- 
plied extensively in the one nucleon sector [20]. In the two-nucleon sector, 
however, the heavy baryon mass formulation of x?T is not applicable due 
to the ill-defined two-nucleon reducible diagrams [7]. Weinberg circumvented 
those problems by deriving chiral power counting rules for time-ordered per- 
turbation theory [7]. 

In this work we follow a different path and work directly with the relativistic 
form of the chiral Lagrangian [1]. This has the advantage that it is possi- 
ble to control covariance at each step of the calculation. We outline how the 



two aforementioned pertinent problems can be overcome within the frame- 
work of the relativistic chiral Lagrangian. We expect our relativistic scheme 
to be equivalent to the more familiar heavy fermion formulation of chiral per- 
turbation theory in the one-nucleon sector. Rather than performing the 1/m 
expansion at the level of the Lagrangian density we suggest to work out this 
expansion explicitly at the level of individual relativistic Feynman diagrams. 
Therewith we avoid the heavy baryon chiral Lagrangian with its known arti- 
facts in the two-nucleon sector. 

It is convenient to consider the chiral Lagrangian (at an intermediate stage) 
to represent a finite cutoff theory [7]. The finite cutoff, A <C m, is required 
to restrict the nucleon virtuality inside a given loop diagram such as to jus- 
tify the 1/m expansion. Obviously a finite cutoff has to be introduced with 
great care not to break any chiral Ward identity. Alternatively one may apply 
dimensional regularization [37,38] . We emphasize that if dimensional regular- 
ization is applied one has to first expand in 1/m and then perform the loop 
integration. The two steps do not commute here. This is equivalent to first 
performing the 1/m expansion and then applying the large cutoff limit in 
a renormalized finite cutoff theory. This procedure suppresses contributions 
from loop momenta larger than the nucleon mass which are the source for 
both the conflict of relativistic diagrams with the counting rule (3) and the 
appearance of multiple powers of the nucleon mass [17]. 

It remains to demonstrate that the relativistic chiral Lagrangian leads to a 
practical scheme in the sense that at a given chiral order Q there are only a 
finite number of diagrams to be considered. Here it is convenient to generalize 
the notation of (3) and denote with di the chiral order of an interaction vertex 
of type %. We will show that by an appropriate regrouping of interaction terms 
in the chiral Lagrangian there are only a finite number of renormalized in- 
teraction vertices of given chiral order. Consider a covariant derivative T>^(x) 
acting on a nucleon field N(x). We need to construct a suitable 'counter inter- 
action' leading to a 'small' renormalized vertex. This crucial rearrangement is 
illustrated at hand of the simplest case when the Lorentz index \x is saturated 
by another covariant derivative acting on the same nucleon field. In this case 
the counter interaction is readily found 



W(x) V^x) N(x) -► (V»(x) V^x) + m 2 ) N(x) ~ Q . (4) 

The 'rearranged' vertex is proportional to the virtuality of the nucleon, (<9 M <9 M + 
m 2 ) N, which is of 'minimal' chiral order Q. It is important to observe that 
this counter interaction is part of the chiral Lagrangian by construction. In 
the more general case where there are 2 k derivatives acting on a nucleon field 
(D^(x) T>n(x)) k N(x) an appropriate rearrangement leads to a renormalized 
structure {T>^{x)V^{x) + m 2 ) k N(x). The total suppression factor of the in- 



teraction vertex di = of + d\ is the sum of the number of derivatives acting 
on pion fields d\ , and the degree of virtuality of the nucleon field d\ = k. 
Note that a similar construction is applicable if covariant derivatives acting 
on different nucleon fields meet. 

We are left to consider the case in which the Lorentz index /i of the covariant 
derivative T>^ is contracted either with the index of a Dirac matrix or the 
index of a derivative acting on a pion field. In the latter case vertices are well 
behaved since they are already suppressed by the pion derivative in d^n. Of 
course, here, one has to give the nucleon derivative the chiral power Q°. The 
'large' nucleon mass is not 'hazardous' if one identifies the natural scale of 
any coupling constant A to be of the same order as the nucleon mass. With 
A ~ m p this seems indeed reasonable. The remaining case in which the nucleon 
derivative couples to the index of a Dirac matrix in a vertex does not generate 
an infinite tower of interaction terms since at a given number of nucleon fields 
rii, the number of available indices provided by Dirac matrices is at most rii. 
The degree of nucleon virtuality can be determined upon inspection case by 
case. 

The above arguments show how to organize the infinite tower of chiral interac- 
tion terms such that there is always only a finite number of interaction terms 
of a given chiral power. The counting rule (3) applies where di is identified 
with the sum of number of derivatives acting on pion fields and the degree of 
nucleon virtuality of the vertex as discussed above. Thus at any given chiral 
order only a finite number of diagrams need to be evaluated. We note that in 
the relativistic scheme it is convenient to use the notion of 'minimal' chiral 
power of a Feynman diagram since a given diagram leads to an infinite tower 
of terms with increasing chiral powers upon expansion. The minimal chiral 
power is given by the counting rule (3). 



2.1 2-Particle reducible diagrams 



In this section we discuss the 2-particle reducible diagrams for which 'stan- 
dard' chiral power counting rules are not applicable [1,7]. In our relativistic 
framework one is therefore bound to consider the Bethe-Salpeter equation for 
the nucleon-nucleon scattering. This is analogous to Weinberg's approach [7] 
which considers the Lippmann-Schwinger equation with a suitable potential 
expanded according to chiral power counting rules. 

The nucleon-nucleon scattering amplitude T, is given in terms of the Bethe- 
Salpeter kernel K 
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where for convenience we use operator notation in the last expression. The 
scattering amplitude T and interaction kernel K of (5) are antisymmetric 
under exchange of the ingoing or outging nucleons. The factor '2' in (5) leads 
to the proper treatment of exchange diagrams. The interaction kernel 




X*X 



+.... (6) 



is the sum of two-particle irreducible diagrams deduced from the relativistic 
chiral Lagrangian. To leading zeroth order it is the sum of local interaction 
terms and the one-pion exchange diagram displayed in (6) with a wavy line 
representing the relativistic pion propagator. The object G in (5) represents 
the relativistic two-nucleon propagator. Let us investigate the Bethe-Salpeter 
equation (5) perturbatively. In the following we introduce the appropriate 
minimal chiral power v for K G K with respect to a selected part of the in- 
teraction kernel K carrying the minimal chiral power vk- The once iterated 
Bethe-Salpeter kernel is 



KOK-r** ^> ™ (7) 

J ( 2vr ) (±w-l) -m 2 + ie(±w + ty -m 2 + ie 

with w 2 = s. In the center of mass frame with w = the Mandelstam variable 
s = 4 m 2 + 4p 2 is given by the relative momentum p of the nucleons. It is 
instructive to study first the s-channel spectral density, p(s), of this contribu- 
tion. It typically receives strength from the pinch singularity generated by the 
two-nucleon propagators in (7). In the center of mass frame one finds: 



p(s) = Z(KGK) ~ -P^ I ^K(t)K(l 

v ' y/m z + p z 4 IT J 4 7T 



Expression (8) can now be used to introduce the minimal chiral power v of 
the s-channel spectral density 

v = 1 + 2 v K (9) 

in terms of the minimal chiral power v^ of the Bethe-Salpeter kernel K given 
by (3). In (9) we exploit that the elastic nucleon-nucleon phase space leads to 
the small momentum p in the spectral density. Note that the angle average in 
(8) does not affect the chiral power. If one assigns K ~ l/(4f 2 ), as suggested 
by the one-pion exchange contribution to the kernel, and identifies m ~ 4:7if n 



the dimensionless expansion parameter p/(4/ w ) arises. This indicates conver- 
gence for small momenta p < 350 MeV. Note, however, that if the kernel K(p) 
in (8) vanishes for large momenta p we expect convergence possibly also for 
p > 350 MeV. This is typically the case if the spectral function p(s) leads to a 
finite loop tftKGK (see (11,41)). We then expect the dimensionless parame- 
ter m 7r /(4 /„.) ~ 0.4 since the falloff in K(l) must be driven by the pion mass. 
It is instructive to compare the chiral power (9) with the 'naive' chiral power, 
2 + 2i>k, implied by the standard counting rule (3). This confirms that the 
counting rule (3) must not be applied to all parts of two-nucleon reducible 
diagrams. We point out that inelastic contributions to p{s), for example from 
the one-pion production cut, can be treated in full analogy to the two-nucleon 
cut. The chiral power of the contribution to the spectral density follows upon 
inspection of the production phase space. In fact here the counting rule (3) 
leads to the correct minimal chiral power (see appendix B.4 and B.5 for explicit 
examples) . 

To have a useful scheme we also need to derive counting rules for the real 
part of the scattering amplitude. Here we invoke causality which relates the 
real part of a given Feynman diagram to its imaginary part by means of a 
dispersion relation 1 



I(s) = $t(KGK) = -V [ ds-^- . (10) 



Am 2 

The dispersion relation (10) need not to be finite and well defined. It may be 
ultraviolet divergent. However, one is free to consider the chiral Lagrangian as 
a finite cutoff theory [7,22] . If the cutoff A is to be interpreted as a bound for 
the maximal virtuality of nucleons one should first remove the nucleon rest 
mass. Therefore we take the cutoff A to restrict the momentum p < A rather 
than the Mandelstam variable s in the dispersion relation (10): 



4(m 2 +A 2 ) _ A 2 

71 J S — S 71 J p A — p z 

Am 2 



1 The iterated Bethe-Salpeter kernel does not necessarily satisfy a s-channel dis- 
persion relation. However, we point out that only the part which provides strength 
for the s-channel spectral density requires special attention and modified power 
counting rules. Obviously a part which does not generate strength for the s-channel 
spectral density has no pinch singularity leading to a s-channel cut so that standard 
counting rules apply. Note that for a given contribution one has to specify whether 
the dispersion relation holds at fixed Mandelstam variable t or u. This technical 
detail, however, does not affect our argument. 



Formally one counts m T /A ~ A/ra ~ Q with m n <C A < m [7]. The natural 
guess for the chiral power of the integral I(p, A) would be the chiral power of 
its imaginary part. However, the integral (11) as it stands does not lead to this 
chiral power yet. Though the separation of scales is manifest in the spectral 
density p(s) the cutoff scale A prevents the conclusion that the integral behaves 
like a small scale to the power 1 + 2 vk- We point out that after an appropriate 
number of subtractions at p 2 = —fi 2 < with fi ~ m n the counting rule (9) can 
nevertheless be applied to the real part of the loop function. The subtraction 
is required to render the dispersion integral finite in the large cutoff limit: 



4»»(p, a;,) = I wV*!£±M (£±4)" . (12) 

71 J p Z — p Z \P + fi J 

The subtraction polynomial can always be absorbed into the contact 4-nucleon 
vertices. Note that the performed subtractions do not necessarily simplify our 
attempt to derive a counting rule since this procedure introduces a new scale, 
fi, the subtraction point. Nevertheless we now can exploit the freedom to 
choose the subtraction point. If we insist on a 'small' subtraction scale fi ~ m n 
of the order of the pion mass one can apply dimensional power counting since 
the integral is finite in the large cutoff limit by construction. 

The above arguments establish the chiral power of a once iterated Bethe- 
Salpeter kernel. The general case of a n-times iterated kernel follows by com- 
plete induction provided that a suitable subtraction mechanism (see section 
2.4) is implemented into the Bethe-Salpeter equation. We arrive at the de- 
sired chiral counting rule for 2-particle reducible diagrams: each pair of nu- 
cleon propagators which exhibit an s-channel unitarity cut receives the chiral 
power —3 in contrast to the 'naive' power —2. Therefore the n-times iterated 
Bethe-Salpeter kernel receives the minimal chiral power 



v k = n + nu K . (13) 

We stress, first, that this counting rule is manifest only if one introduces a 
'small' subtraction scale /i ~ m^ at the level of the Bethe-Salpeter equation 
and second, that the subtraction scale is independent of the intrinsic cutoff. 
We will refer to this counting rule as the 'L '-counting rule since for the one- 
pion exchange interaction, which contributes at order zero to the interaction 
kernel, the chiral power of the n-th iterated interaction is simply given by the 
number of loops. In appendix B the reader may find explicit examples which 
confirm our L-counting rule (13). 

We summarize what we derived. The non-polynomial part of the once iterated 
Bethe-Salpeter kernel has an excess of one chiral power as compared to the sum 



of the chiral powers of the iterated Bethe-Salpeter kernels (see (9)). We point 
out that this result, though obtained by means of dispersion relations, holds 
also if a different technique for the evaluation of the loop functions is applied. 
Since the chiral power of the Bethe-Salpeter kernel starts at v — our counting 
rules seem to suggest that the Bethe-Salpeter kernel has to be iterated only a 
finite number of times at a given chiral order v for the scattering amplitude. 
Obviously this can not be correct since there are well known non-perturbative 
effects present in the scattering amplitude like for example the deuteron bound 
state. This puzzle can be resolved by realizing that the implicit assumption of 
any chiral power counting scheme is the naturalness of all coupling constants 
carrying non trivial dimensions. The characteristic scale is assumed to be the 
mass of the lightest degree of freedom not explicitly included in the effective 
chiral Lagrangian. 

Therefore the crucial question to ask is: does the suggested subtraction scheme 
for the Bethe-Salpeter equation affect the naturalness of the J^-nucleon coupling 
strengths ? 



2.2 Mutation of natural scale to anomalous scale 



In this section we address the question to what extent a subtraction scheme 
affects the naturalness of the bare coupling function. The effect of our sub- 
traction scheme is illustrated at hand of s-wave scattering in the isospin one 
channel. We consider here a theory without explicit pions defined by an infinite 
tower of separable interactions in the 1 S'o-channel 2 : 



-I oo 

£an(x) = - J2 9 {2n) d n (N 75 C tt 2 N f ) 8 n (TV* r 2 r C' 1 ^ N) (14) 

4 n=0 

with the charge conjugation matrix C = 27072 and the isospin Pauli matri- 
ces Tj. The index contractions of the derivatives d^ in (14) are such that in 
momentum space the energy dependent coupling g(s) arises: 



00 00 

</(*) = E 9 (2n) s n = E 92n{s-Am 2 ) n (15) 

n=0 n=0 



2 Note that the terms displayed in (14) are the only ones relevant for s-wave scatter- 
ing in the spin singlet channel. Any further terms necessarily are on-shell equivalent 
to terms already included and lead to relativistic nucleonic tadpole diagrams upon 
unitary iteration in the Bethe-Salpeter equation. Such terms can be dropped in a 
consistent treatment (see appendix B.l and (46)). Note that the situation is different 
in a non-relativistic framework where one encounters a richer variety of interaction 
terms [26]. 



10 



with p 2 = s/4 — m 2 . In (15) we introduce renormalized couplings Q2 n accord- 
ing to the mandatory regrouping of interaction terms discussed in section 2. 
Note that in the derivation of Feynman rules from (14) the troublesome time- 
derivatives can be treated naively since the form of the final effective Feynman 
vertex is dictated by covariance. At tree level the coupling go can be interpreted 
in terms of the s-wave scattering lengths a( 1 S'o) with 2mgo = — 4 7ra( 1 S'o). Re- 
call that one may apply a phenomenological saturation model for the coupling 
function with g(s) receiving contributions from heavy meson (cr,p,uj, ... ) ex- 
changes. Such a model provides a bridge to the traditional meson-exchange 
picture of nucleon-nucleon scattering. It furthermore motivates the crucial as- 
sumption of any chiral expansion, namely, that all couplings are natural: in 
particular we expect g n+ 2 ~ <?n/A 2 at . We identify here the natural scale A nai . 
with one half of the mass of the lightest meson not included explicitly in the 
theory. The factor 1/2 follows since the t-channel meson-exchange diagram in 
a given partial wave suggests a convergence radius of half the meson mass for 
an expansion around p 2 = 0. In a theory with pions one identifies A nat . = m p /2 
and in a theory where pions are integrated out A nat , = m n /2. It is worth em- 
phasizing that the leading term requires special attention since g ^ 1/A^ at . 
For example the vector meson u and t channel exchange diagrams predict 

9o 9 2 P NN/( 8m2 P ) -9l N N/( 8m l) ~ ~ 1 /( 8 f 2 ) with a typical meson-nucleon 

coupling constant g^NN ~ 9 p nn ~ m p /(y/2 f w ) estimated with the KSFR re- 
lation. Therefore one expects go = A 2 go with the dimensionless parameter go 
of order unity. 

We proceed and solve the Bethe-Salpeter equation for the nucleon-nucleon 
scattering amplitude T, with the tree-level Bethe-Salpeter kernel of (14). The 
equation is separable and therefore the solution can be written in terms of the 
divergent one- loop integral, J(s,A), 



J(s, A) = -i tr J ^- CT'ts S f (l + \w) T 5 CS t F (-l + \w) . (16) 

Here Sf(p) = l/(rf> — m + ie) is the relativistic nucleon propagator and 
w 2 = 4(m 2 + p 2 ). The loop integral may be evaluated by writing a disper- 
sion integral 3 



3 We are aware that the introduction of a finite cutoff in Feynman diagrams requires 
special care in order not to break the pertinent symmetries. For more details we 
refer to the appendix. 
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J(S , A) = I/4^ L (17) 

it J p 2 — p 2 — ie 



where we restrict the relative momentum of the center of mass system by the 
finite cutoff parameter A. The well defined spectral density is: 



2tt \ 2 m 2 

We expand the loop integral J(s, A) non relativistically in powers of 1/m by 
expanding its imaginary part (18). Here we encounter a simple example how 
to systematically expand a Feynman diagram given by the relativistic chiral 
Lagrangian. Consider for example the loop function at zero momentum 



A 





J(4m 2 ,A) = — - A \Jm 2 + A 2 + m 2 arcsinh — . (19) 

2ir 2 \ \mj J 

A naive evaluation of the diagram with A ^> m would generate a troublesome 
multiple power of the nucleon mass with J ~ m 2 log A counteracting any chiral 
counting rule. Equations (17) and (19) exemplify the necessity of the intrinsic 
cutoff A to be smaller than the nucleon mass such that the 1/m expansion 
is justified. Expanding (19) in A/m reduces the leading power of m. This is 
also seen if the real part of the loop J(s, A) is evaluated with the expanded 
spectral function (18) by means of the dispersion integral (17) : 



m I A 

MJ(s,A) = — A 1 + 




0(p 4 ) . (20) 

It is worth pointing out that the subleading terms in the 1/m expansion 
proportional to (A/m) n with n > 1 can be evaluated with A = m. This follows 
since the A/m expansion is controlled typically by square root singularities 
y/m' 2 + p 2 for which even at A/m = 1 the non-relativistic expansion (18) is 
converging rapidly. This is nicely confirmed by the suppression factors 1/6 
and 1/12 for the constant term and p 2 -term in (20) respectively. Therefore, 
effectively we are left with only linear divergencies in A. In fact it is legitimate 
to drop all terms (A/m)™ with n > 1 since at given order p 2k their effect 
can be compensated by a small change of A. This is particularly convenient 
since then nucleonic tadpole diagrams can be dropped systematically. With 
this simplified regularization scheme we avoid the problems one encounters in 
cutoff-regularized effective field theory. The ill-defined multiple powers of the 
cutoff are eliminated and need not to be canceled by the careful construction 
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of the cutoff dependence of path integral measure [27]. For more details on 
our regularization procedure including or more rigorous definition we refer the 
reader to appendix B. Here we point at a further important observation to be 
made. The presence of terms like (p/A) 2k necessarily leads to the identification 
A ~ A nat .. Only then such terms can be absorbed into contact interaction 
terms of the chiral Lagrangian. 

To leading order in the non-relativistic expansion the one-loop integral (17) 
can be rendered finite by one subtraction at s = 4 m 2 — li 2 



A 2 

Ms, A; vl) = J(s, A) - J {A m 2 - ^ A) = I / _f * p ^\ t±^ (21) 

it J p 2 — p 2 — ie p 2 + ii l 

where the subtraction constant J {Am 2 — li 2 ,A) can be absorbed into the cou- 
pling function g(s) = g(s, A). Note that subleading terms in the 1/ra expansion 
require further subtractions. The subtracted loop (21) is not identified with 
the renormalized loop. It is more convenient to deal with a manifest cutoff 
independent renormalized loop function. We introduce Jr(j), li) = Js(s, oo;/x) 
and accordingly the renormalized coupling gii(p 2 ,[i) with 

T[1Sa] {P) = g-i(s,A)-J(s,A) = g R \p 2 ,n)-J R (p,n) ' (22) 



All dependence on the way in which the cutoff A was introduced is absorbed 
into the bare coupling function g(s,A). Consequently the choice of the cutoff 
scheme does not have any effect as long as the scheme does not break the 
pertinent symmetries. The leading order renormalized loop reads 



oc 



m r dp p p ' + /j, m 

2n 2 J p 2 — p 2 — ie p 2 + fx 2 2n' 



Jr(p,») = -Z-2 / ^2_^_,-^2 I l|2 =7^(v + ip) (23) 





where we dropped all 1/m correction terms in (23). The leading order renor- 
malized coupling function is: 

-1/ 2 -1/ a\ m f dp 2 p m f dp 2 p p 2 + fi 2 



9r(P,H)=9 ( g » A ) ~ 7T^ / ^ 2 I , 2 + 



2n 2 J p 2 + [x 2 2ir 2 J p 2 — p 2 p 2 + ji 2 

g-\s, a)-™U-1 im ) + ™?(i + o( p 2 /a 2 )) . (24) 



IX 2 \ Z / 77 



A 



It is important to observe that for p < A the contribution form the loop func- 
tion to the renormalized coupling function can be expanded in powers of the 
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small momentum p. Our result (24) demonstrates that all cutoff dependence 
can safely be absorbed in the bare coupling g(s,A) provided that it is legit- 
imate to identify the typical cutoff A with the natural scale of the bare cou- 
pling function A nat .. Equation (24) exhibits an interesting effect. Starting with 
a natural bare coupling function g(s,A) it potentially generates an anoma- 
lously large renormalized coupling gn(p 2 = 0, //), however, only if the subtrac- 
tion scale \x is taken to be sufficiently small as required by the L-counting 
rule. Then an attractive natural coupling g~ l (4m 2 , A) > potentially cancels 
the mA term in (24) predicting an anomalously large renormalized coupling 
gn(0,n) ~ 27r/(m/z). We stress that this is a familiar phenomenon since it 
is long understood that an attractive system may dynamically generate 'new' 
scales which are anomalously small (like for example the binding energy of 
shallow bound state). What is intriguing here is the transparent and simple 
mechanism provided by (24) to naturally generate such a 'small' scale [1,28]. If 
the cancellation in (24) is manifest for a physical system the renormalized cou- 
pling gn(0, //) must be assigned an anomalous chiral power, e.g. the 'minimal' 
chiral power —1 at given subtraction scale with \x ~ m 7r . 

In (24) we refrained from expanding the inverse bare coupling function g~ 1 (s, A). 
We point out that it is necessary to discriminate two scenarios. In scenario I) 
the convergence radius of g~ 1 (s, A) is set by the appropriate natural scale A„ ai 
leading to a well behaved Taylor series in p 2 . In scenario II) the convergence 
radius of g~ x (s, A) is set by a zero of the bare coupling function at p 2 = fj, 2 with 
/j, a < A„ at . In this case it may be advantageous to expand the renormalized 
coupling around p = \x a . If one insists on expanding around p = the theory 
breaks down at scale p < fi a < A„ at . In fact expanding around p = /i CT further 
increases the convergence radius with A„ at . = A„j 4 + /i CT simply because the 
distance of the expansion point p = /i CT to the closest subthreshold singularity 
at p = ± i A ( n at. is increased. Note that here we implicitly assume that /v+A„ at 
is smaller than the one-particle production cut with 2 /i^ + /i CT A„ at < m A„ at . 
For example, in an effective theory which treats pions explicitly and where 
the p-meson is the lightest relevant degree-of-freedom that is integrated out, 
the optimal expansion point p ~ —m p /A + [mm p /A + vn 2 jVo) x l 2 ~ 275 MeV 
leads to the maximal convergence radius or equivalently a natural scale of 
A nat . ~ 660 MeV. 

We investigate the scattering amplitude T|is ](p) of (22) in the two scenarios 
to leading order in the 1/m expansion 



m ■ " ' 



T [ ( 4]G°) = {gh~%) + 9 { rj 1] p 2 + ■ ■ ■ - ^ {» + ip 

\P Av 
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where we implicitly assume that the bare coupling shows a simple zero with 
g(/J%.) = but g'(iJ%) 7^ 0. In principal a zero with higher degeneracy can not 
be excluded though it is less likely for a natural coupling function. Multiple 
zeros of g(p 2 ,A) are difficult to reconcile with the naturalness assumption. 
We emphasize that an effective field theory formulated with an explicit cutoff 
may give rise to a similar behavior since here the inverse bare coupling is 
kept and not expanded. The zero at p 2 = fi 2 , however, is then necessarily a 
consequence of strongly correlated and fine-tuned parameters resulting from 
a fit to data. In our scheme (25) we avoid a strong correlation of parameters 
since the vanishing of the phase shift is already build in and controlled directly 
by /i CT . There is a further important difference to our scheme: at order Q 2 the 
inverse scattering amplitude of scenario I) requires two parameters whereas 
scenario II) leads to four parameters. In other words scenario II) predicts the 
relevance of more operators at a given order. 

It is instructive to explicitly work out the implications of the naturalness 
assumption and discuss the low energy structure of T^s ]{p) at momenta p <C 
\x a in both scenarios. It leads to coupling functions of the form: 



9^ M = ~ (Ail -?**)+ 9%r 4 fl 

2 



9 \ 1HJLL. c\ 

7T V 2 






where the dimensionless quantities g^j'f , g^n and g R 'j are of order 
unity. We point out that g R ' n ~ 4 f 2 and g R 'j (/x) ~ 4 f 2 holds only for 
fi ~ A„j t and /j, ~ h- na t. respectively. The subleading coefficients g R 'f and 
9rii d° no ^ depend on the subtraction scale // for n > 2 to leading order in 
the 1/ra expansion. Therefore, at fj, — the coefficient —g RI } /^a+dRii (0) 
and also g RI (0) may be anomalously small leading to an unnatural large 
scattering length 2ir a = — mTji5 ](0). The required cancellation is achieved 

with an attractive coupling g R j ^(0) ~ A„j t ./(7r/ ff ) in scenario I) and with 

g RI f — —fiff/i^fir^iiat.) i n scenario II) where we identified 471/7,- ~ m for 
simplicity. 

An anomalously large scattering length is evidence for the presence of non- 
perturbative structures at small momenta. If the precise singular structure is 
known one may proceed and break up and further expand the denominator of 
(25) without loosing accuracy in the expansion. Note, however, that such an 
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expansion will necessarily be in conflict with unitarity. In scenario I) a system- 
atic expansion of the denominator is induced by identifying g R 'j (m„-) ~ Q 

and g R 'j (m„-)p 2 ~ Q 2 with the typical small scale Q. This implies that the 
scattering amplitude T[is ](p) shows either a shallow bound state or a pseudo- 
bound state [1,15]. The perturbative inclusion of pion dynamics according to 
the L-counting rule is straight forward. If the scattering amplitude Tng i(p) 
exhibits a more complicated low-energy structure the scenario I) must be re- 
jected as inefficient with an unnaturally small convergence radius. In scenario 
II) the relative importance of the two terms g R jj and g R ' IT (/J,) depends 
crucially on the subtraction scale /x. For small subtraction scale \x ~ m n , re- 
quired for a manifest realization of the L-counting rule (13), both terms are 
equally important at p ~ 0. In fact the low-energy characteristics of scenario 
II) is much richer than that of scenario I). At 'leading' orders the inverse 
scattering amplitude may have zeros at three different complex momenta cor- 
responding to the roots of a cubic equation. It is unclear, in particular when 
the pion dynamics is included, in which way to expand the denominator with- 
out loosing the homogenous convergence property of (25). In scenario II) it is 
therefore advantageous to keep unitarity exactly and refrain from any further 
expansion beyond (25). 



2.3 S-wave scattering in the spin singlet channel 



Before evaluating Tng i (p) to leading orders in the chiral expansion we recall its 
empirical behavior. A comparison of the structure of the empirical scattering 
amplitude with that of the amplitude in the two scenarios (25) provides strong 
indications to which scheme to choose. 

The empirical scattering amplitude shows two striking phenomena: a pseudo- 
bound state close to threshold and a vanishing phase shift or equivalently a 
zero of the scattering amplitude Tng -\ = at p ~ 355 MeV. At very low 
energies with E\ ah < 5 MeV the partial wave amplitude T[i 5o ](p) is dominated 
by a pseudo-bound state pole 



T [ l S ]iP) = -7= 




»tih>+0[£\) (27) 



p + iz ' \m 



where we made the pseudo-bound state pole contribution explicit and ex- 
panded the remainder in powers of p. The pole position z, its residuum, w, 
and the leading background parameters b e $ are predicted accurately by the 
effective range theory 
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1 1 

p cot 5 = — + - rp 2 + O (p 4 ) (28) 

in terms of the s-wave np-scattering length a ~ —23.75 fm and effective range 
parameter r ~ 2.75 fm. One finds 



r t 1 — w 



w- 1 = \/l-2-~0.90- 1 , z = ~0.04fm" 1 , 

V a wr 

,(q) rw 2 w 1 — w 

oS'^ + T^-iTT^- 1 - 18671 ' 

$ = »_„.= » (IZH^IW (29) 

where we included also the prediction of the effective range theory for the 
next-to-leading moment in (27) for completeness. Note that for example the 
effective scattering length b c{ { is substantially suppressed as compared with the 
empirical scattering length a ~ —23.75 fm 4 . We point out that the effective 
range formula (28) implies a second pole in the scattering amplitude at p ~ 
% 0.77/fm: 



T ?sM 



cflf. range 



W 




IW 



p — I z 



- z = \±^ z= l + 0[-\ . (30) 



From a theoretical point of view it is not obvious that the second pole at 
p = iz carries direct physical significance since it may 'move' as more terms 
are included in the expansion of the inverse i^-matrix. In other words it could 
be an effective singularity which accounts for non-polynomial contributions 
from the exchange of pions. Note, however, that from scattering theory one 
expects that the analyticity domain of the s-wave amplitude is bounded due to 
the one-pion exchange contribution at Qp < m n /2 and m^/2 < z numerically. 

We now discriminate two alternative scenarios. Suppose first that the strength 
of the pole at p = i z effectively accounts for perturbative pionic cuts. Then 
one may expand the amplitude as indicated in (27). The chiral theory is ex- 
pected to replace the polynomial remainder in (27) by smooth and calculable 
analytic functions which exhibit the proper pionic cuts. A natural starting 
point for a theoretical description is scenario I) of (25) with the denomina- 
tor in T,[g Ap) expanded according to the L-counting rule. This leads to a 
perturbative treatment of pions. 



4 It will turn out to be more convenient to work with the parameters w and z. 
The scattering length a and effective range are recovered as follows: r = ^j- and 

n - 2w . 
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Alternatively, suppose that the pole at p = iz has a physical significance 
in the sense that the exact solution of the problem shows an anomalously 
large subthreshold amplitude in the vicinity of p ~ % z. Then the expansion 
(27) turns useless at rather small momenta p < z. Therefore an effective field 
theory description must incorporate that singularity properly. A promising 
starting point is now scenario II) of (25). We expand the inverse of TL, 
and identify /x CT with the empirical zero of the scattering phase. In this case 
we expect a significantly larger natural scale A; iat { > p, a ~ 350 MeV. In the 
next section we show in detail how the pion dynamics can be incorporated 
perturbatively according to the L-counting rule. 

In this paper we will argue that the second scenario is in fact realized in 
nature. The success of the effective range theory, which describes the 1 Sq 
phase accurately up to momenta p ~ 150 MeV, together with the L-counting 
rule already leads to the second scenario. Since the effective range theory can 
be viewed as an effective low energy theory where the pion degrees of freedom 
are integrated out, the second pole at p ~ % 150 MeV points either at the 
importance of unitary two-pion exchanges with a branch point at p — % m n 
or at the physical significance of the second pole. On the other hand the 
chiral L-counting rule predicts that the pionic cut at p — % m n is weak. From 
the fact that the effective range theory works well up to surprisingly large 
momenta p ~ z ~ m n and the effective range theory prediction of a strong 
second pole at p = iz with residuum w ~ 1 (see (30)) we then conclude 
that we have choose the second scenario. This leads to an anomalously large 
subthreshold amplitude in the vicinity of p ~ iz One arrives at a similar 
conclusion upon inspection of the typical size of the effective range parameter 
in both scenarios. Scenario II) yields a natural interpretation of the rather large 
empirical effective range parameter \jz ~ r ~ 2.75 fm also in a model without 
pions. We return to this point in section 4 where we show that the empirical 
scattering amplitude can be rather nicely described by T\ s \ip) parameterized 
in terms of the four leading parameters in (25). 



2.4 Solving the Bethe-Salpeter equation 



In the previous section we derived a chiral power counting rule for appropri- 
ately subtracted 2-particle reducible diagrams generated by the once iterated 
Bethe-Salpeter kernel. In order to show that the counting rule (13) holds for 
more than one iteration of the Bethe-Salpeter kernel we now elaborate on how 
the subtraction scheme is implemented in the Bethe-Salpeter equation 



T = 2K + KGT . (31) 



The interaction kernel K is the sum of polynomial contributions from 2- 
nucleon contact interaction vertices, K a , which represents the exchange of 
heavy meson (a,p,u,... ) and the term K n , which represents 2-particle irre- 
ducible pion exchange diagrams: 




X*X 



+ ... = K <T + K ir . (32) 



Obviously the above separation of the interaction kernel K is not unique since 
the irreducible multi-pion loops may also generate polynomial contact inter- 
action terms. We will elaborate on this issue in detail below. Let us introduce 
the amplitude T n 



T„ = 2(l-K«G) X K, (33) 

which is the solution of the Bethe-Salpeter equation for K a = and the 
auxiliary scattering amplitude T a : 






2 A'„ 



k - ] -^KCS 



2 ^ _ 

-i 



= 2(\-K a G-\K a GT^G) K a . (34) 

The formal solution of the Bethe-Salpeter equation is easily expressed in terms 
of these amplitudes T a and T n 



m 






T) - (T.) + - j^J^ + 2&USQ 



l + ±T,G)Tjl + \GT,) . (35) 



Next we introduce a subtraction scheme for the Bethe-Salpeter equation. The 
subtracted kernels K n ^$ and K a g are defined in terms of an arbitrary and not 
yet specified subtraction interaction K$ 

K K)S = K n + K s , 

K U , S = K U -K S . (36) 
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According to our construction the full interaction kernel K and consequently 
also the full scattering amplitude are manifestly independent of the subtraction 
interaction K$ 



K = K 7r + K a = K VtS + K^ s . (37) 

The subtracted interaction kernels induce the subtracted scattering ampli- 
tudes T n> s an d T u>s with: 



T^ s [K s }=2(l-K^ s Gy 1 K^ s , 

T <TtS [K 3 ]=2(l - K^ s G-\K^ s GT^s[Ks\ G)~' K a , s . (3* 

The full scattering amplitude remains formally invariant 



T = T n , s [K s \ + (l + \ T VtS [K s ] G) T a , s [K s ] (l + \GT^ S [K S \) . (39) 



One is now free to exploit the choice of the subtraction interaction Kg- Differ- 
ent choices for Kg may increase or decrease the average interaction strength in 
K n ^s- Therefore one would expect that the non-perturbative part of the scat- 
tering amplitude like for example the deuteron bound state may sit either in 
T^ s or T n> g depending on the choice of Kg. Suppose now that non-perturbative 
effects sit in T^g then one may conclude that the amplitude T^g is perturba- 
tive and can be expanded in powers of K n ^g. We point out that in this case 
also the amplitude T a> g need not to be evaluated exactly. The perturbative 
character of the subtracted pion iteration permits the expansion of the inverse 
amplitude T ff s m powers of the kernel K„ t g. 



T a , s = 2 (l - K ff , s G - K ff , s G K v , s G + ...) V^ 

= 2 (k-1 -G-GK V>S G+ ...r 1 . (40) 



Here we refrain form breaking up and further expanding the denominator 
in (40) since that requires further assumptions on the structures of non- 
perturbative effects in T aj g. 

It is useful to outline our renormalization procedure. At given chiral order n, 
construct the local kernel Kg = Kg perturbatively such that 



T KtS [K s ] = finite + O (Q" +1 ) & T^ S [K S ] G = finite + O (Q n+2 ) (41) 
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holds. Note that T^g G can be defined via T n> s G K a s if K^s is separable. We 
emphasize that the construction (41) differs from conventional renormalization 
procedures insofar that the object T w ^ s G is rendered finite by 'counter loops' 
rather then counter terms. In particular T^g and T^g G do not depend on 
any renormalization scale. As an immediate consequence note that for T n> s 
and also the vertex loop functions T n> s G the L-counting rule apply without a 
further subtraction procedure (see (12)). The only remaining divergencies sit 
in T a s and can be absorbed into the local bare coupling K a : 



K;} -G- GT KtS [K 3 ] G = finite + O [Q n+A ) . (42) 

Note that the divergencies in K a s G T^g G K a ^s are 'weak' since they are 
already tamed to a large extend by the renormalization condition (41). Here a 
manifest realization of the L-counting rule requires a subtraction scale \x ~ m w . 

It is not obvious that the renormalization conditions (41) and (42) can be met. 
In fact our scheme requires that the Bethe-Salpeter interaction kernel K is 'on- 
shell irreducible' if we insist on a separable counter interaction K$- Basically 
we need a mechanism which absorbs systematically all divergencies not directly 
linked to a s-channel unitarity cut into K a at the level of the interaction kernel. 
We decompose the kernel K of (31) into an 'on-shell irreducible' part K and 
'on-shell reducible' terms Kl, Kr and Klr which vanish if evaluated with 
on-shell kinematics either in the incoming or outgoing channel 



K = K + K L + K R + K LR . (43) 

Here Kr (Kl) vanishes at on-shell kinematics for the incoming (outgoing) 
nucleons after multiplication with the incoming (outgoing) nucleon spinors 
by construction. The term Klr vanishes if evaluated with either incoming 
or outgoing on-shell kinematics. For example the s-channel spectral density of 
Kr G K or K G Kl is zero. A term like Kr G would be renormalized according 
to (41) by an appropriate counter loop K$G with a non-separable K$- This 
would destroy the simple structure of T^ s . It is therefore advantageous to 
eliminate all on-shell reducible terms by introducing the on-shell equivalent 
effective interaction kernel K with 



T = 2K + KGT = - —-K 

I- KG on-shell 

K = K + K r GK + KGK l + KrGK l + KGK lr GK + - . (44) 

By construction the effective interaction kernel K is on-shell irreducible and 
one therefore expects a non-zero s-channel spectral density of K G K. We point 
out that it is sufficient to work out the matching in (44) perturbatively since 
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chiral power counting rules apply for K 5 . In full analogy to (37) the effective 
kernel K = K n + K a is written as a sum of non-local contributions K n and 
local counter terms K a . Divergent terms in K can now be absorbed into the 
local interaction vertex K a of (32) at the level of the interaction kernel. 

We need to address a further aspect of our scheme: how does one renormalize 
divergencies that appear in sub leading orders in the 1/m expansion. To this 
end it is convenient to introduce a partly renormalized two-nucleon propagator 
G: 



4(m 2 +A 2 ) _ ,— — 

G(l;w)= f —\— 2 ^— AS { +\l,w;s) ® AS { f\l,w;s) , 

J 7i \ w z s — w z — t e 

Am 2 

AS { ^\l,w;s) = iTrS(l 2 ±l-w + \s-m 2 ) (±/+^ + m) (45) 

where we keep the cutoff parameter A at this intermediate stage. Replacing 
the two-nucleon propagator G by the renormalized propagator G leads to a 
renormalized interaction kernel K 



f = 2k + KGf = 2K + KGf, K = ^ —K. (46) 

l-K{G-G) 

It is crucial to observe that by construction K is free of s-channel unitarity 
cuts. Therefore K can be evaluated in perturbation theory applying chiral 
power counting rules. The construction (46) defines a systematic prescription 
for absorbing ill-defined power divergent terms (A/m) n with n > 1 induced 
by the 1/m expansion into K a (see (20)). Note that for a local vertex (14) we 



would find K = K if we replaced Js/w 2 — > 1 in (45). 

In the following we assume that the kernel K in (31) is on-shell irreducible and 
1/m renormalized according to (44,46). We emphasize that the constructions 
(44) and (46) are needed to systematically absorb a large class of divergencies 
into the local interaction kernel K a . Only after the effective interaction kernel 
has been properly identified, can one apply the renormalization conditions 
(41,42). Note that in actual calculations it is not mandatory to construct the 
effective interaction explicitly. Since the pion-dynamics at given chiral order 
leads to a finite number of Feynman diagrams to be considered it is sufficient 
to evaluate those diagrams in terms of the original bare interaction kernel 



5 Note that the introduction of the on-shell irreducible kernel via (44) demonstrates 
that the infinite tower of local s-wave interaction terms in (14) is effectively complete. 
Further s-wave interaction terms simply renormalize the coupling function of (14) 
according to (44). 
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K. The construction (44) and(46) then leads to a prescription for moving 
divergencies, which are not absorbed with (41,42), into K a (see appendix B). 

We turn to a problem with the proposed expansion: as it stands it violates 
unitarity. Here we anticipate an important result to be discussed in detail 
below, namely that such an expansion can be reconciled with unitarity if 
the subtraction interaction Ks is generalized to allow for more complicated 
separable structures. Rather than evaluating for example T^ t s[Ks] strictly at 
Ks = Ks as suggested by the renormalization condition (41) we consider 
Tjt^sIKs] at Ks ^ Ks- It is crucial not to spoil the renormalization procedure 
(41,42). Therefore the non-trivial part of Ks is permitted only in the finite 
parts. More specifically we split for example T w ^s[Ks] m a finite part T* l g'[Ks] 
and a divergent part T^ l g-[Ks] and expand in Ks around Ks'- 



T w ,s[p 2 ; Ks] = Tl m s \V\ tf\ K 8 ] + T«f [p\ A* 2 ; K s ] 

oo 

= E (?tL(p 2 , f) + Tl% n {p\ i* 2 )) {Ks -Ks)\ 

n=0 

Tl% n {-U\j?) = if Ti;L(p 2 ^ 2 ) ± (47) 

where we introduce a necessary splitting point \x. In (47) the local kernel Ks 
is fixed according to the renormalization condition (41). We implicitly assume 
a projection onto a given partial wave. Thus the scattering amplitude is a 
function of only the relative momentum p with s = 4 (m 2 + p 2 ). Note for 
example that by construction T^' div (p 2 , Jx 2 ) = and therefore T^j in (p 2 , p 2 ) is 
independent of the splitting point p,. Moreover in (47) we assume that the non- 
trivial part of the kernel Ks is separable since only then the multiple powers 
of (Ks — Ks) n are defined properly. The separation of the objects T n> s G and 
GT„ } s m to finite and divergent parts is performed in complete analogy with 
(47). The remaining objects GT n ^G and K a s require a separate treatment 
since they involve further divergencies which are not taken care of by Ks'- 



K- s -G = GT„ :S [Ks]G 

oo 

+ E (rJX,! V, f) + t^( p 2 , f)) (k s - k 



T^XK-fi 2 , f) = if T^d?, f) ± . (48) 

Again by construction T^' d ~ v '(p 2 ,JX 2 ) = and therefore T^'j^ n '(p 2 ,jj, 2 ) is in- 
dependent on the splitting point Jx. Note that the splitting procedure in (47) 
and (48) must necessarily be correlated in order to ensure that the full scat- 
tering amplitude T is independent of Ks and thus on any power (Ks — Ks) n 
as it should be according to our construction (36,37). The independence of 
T on Ks — Ks serves as a powerful argument that all divergent terms like 
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T^ div (p 2 , n 2 ) [K s — K s ) can be dropped. Note that any truncation of T^, 



and T CTj 5 of the form (40) necessarily leads to a residual dependence of T on 
Ks- This dependence of the approximate T on Kg — Ks can then be exploited 
to restore unitarity. 

We proceed and present an argument that the splitting point p 2 must in fact 
be identified with the chiral limit value of a potential singularity present in 
T at p 2 = —z 2 = —fi 2 in a given partial wave 6 . Suppose we drop the non- 
trivial terms (Ks — Ks) n for a moment since in the full expression for T they 
cancels anyway. Then we expect the singularity to sit in T„j in (— Zq) = since 
T^ \ in (p 2 ) is perturbative and independent on p, by construction. Therefore, if 

we now include the terms (Ks — Ks) n , we must insist on T^j in (—z 2 ,j2 2 ) = 
also for n > 0. If we do not insist on this condition the singularity at 
p 2 = —z 2 must sit in T^ l in which contradicts our construction. Here we 
implicitly assume that the non-trivial part in Ks is also perturbative. This 
will be confirmed later when this part is constructed explicitly. We now arrive 
at our result fi 2 = z 2 , since from (47) it follows that T^ iv '~ (p 2 ,j2 2 ) ^ holds. 

We emphasize that the local part of Ks is fixed by the renormalization con- 
dition (41) whereas the non-trivial terms (Ks — Ks) n , which are included 
exclusively in the finite parts, are determined by the unitarity constraint. In 
our discussion of the separation of finite parts from infinite parts, we sup- 
pressed the chiral expansion of e.g. Ks and T ffjl g. It is implicitly understood 
that (47,48) apply order by order in the chiral expansion. 

On a rather formal level the underlying mechanism leading to exact unitarity is 
quickly revealed. At given chiral order n construct the non-trivial but separable 
part of Ks such that 



(T^s) = + O (Q n+1 ) (49) 

holds, where the brackets refers to projection onto a given partial wave with 
on-shell kinematics. Insisting on a separable construction of Ks leads to the 
crucial identities: 



(1 + \T^ S G) (1 + jGT„ tS ) + O (Q"+ 2 ) 
(K-],-G-\GT^ s G) + 0(Q^) ' 

= ^(l+ l -T^ s G)(l+ l -GT^ s ) + O{Q n+2 ) . (50) 



6 In the spin triplet channel for example we would identify zq = J 2 meo — e 2 D with 



-D 

the nucleon mass m and the deuteron binding energy cd evaluated in the chiral limit 
with 77V = 0. 
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The unitarity condition 



QT- 1 = ~$fG (51) 

applied to (39) reduces to: 



3 (G + \ GT^ S G) = (3 (G)) (ft (1 + \ T v , s G)(l + ± GT^ S )) . (52) 

The unitarity condition (52) is indeed fulfilled if expanded according to chi- 
ral power counting rules. This is basically Cutkosky's cutting rule [23]. The 
explicit technical realization of this mechanism will be presented in section 
3.4. 

According to the chiral counting rules one expects that any non-perturbative 
structure must sit in the denominator of (50). Thus any attempt to break 
up and expand further this denominator is subtle if the precise singularity 
structure is unknown. We suggest a scheme in which the numerator and de- 
nominator of (50) are expanded according to the counting rule (13) separately. 
This expansion has the great advantage that it respects unitarity exactly order 
by order in the chiral expansion. Note that if T n< s = the full solution of the 
Bethe-Salpeter equation is represented by T^s- In this sense one may say that 
unitarity leads to an optimized chiral expansion of the scattering amplitude. 
We stress that a perturbative K n> s kernel is predicted by chiral power count- 
ing rules provided that the subtraction interaction Kg is constructed properly. 
In particular all divergencies must be absorbed in T^g and T n> g G. A positive 
feature of this scheme is its practicality: the object T^g can be evaluated per- 
turbatively and therefore, due to the separable nature of Kg, the scattering 
amplitude can be evaluated analytically. 



3 The nucleon-nucleon scattering amplitude 



In this section we derive explicit expressions for the leading chiral orders of 
the nucleon-nucleon scattering amplitude constructed according to the chiral 
counting scheme of section 2. We will present analytic results for the required 
loop functions and confirm the usefulness of our L-counting rule. 

We parameterize the on-shell scattering amplitude, T(t,u), in terms of five 
isospin zero amplitudes, a] (t,u), and five isospin one amplitudes, a\ (t,u), 

5 5 

T(t, u) = VoYl aS 0) ^ u ) T l + V 1 Y1 4 1} (t, u) T % (53) 

i=0 i=0 
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defined with respect to the Dirac tensors 



7i = 


(75 C-'xC 75) , 


T 2 = 


( 7a C- 1 xC 7 a ) , 


T 3 = 


(C _1 x C) , 


T 4 = 


(7.C- 1 x C 7 a ) - \ (o^C- 1 x C/) , 


T,= 


(75 C- 1 x C75) + (757a C" 1 x C l5l a ) 



(54) 

with the charge conjugation matrix C = 27072 and the isospin projection 
operators 

V = \ (r 2 x r 2 ) , 

7>i = \ (rr 2 x r 2 r) . (55) 

For the precise definition of the tensor product x we refer to appendix A. 
The invariant amplitudes, a\ (t,u), depend on the Mandelstam variables t = 
(pi — Pi) and u = (pi — p' 2 ) 2 where pi j2 are the 4-momenta of the incoming 
nucleons and p[ 2 those of the outgoing nucleons. The Dirac tensors Tj are 
constructed such that they are invariant modulo a sign upon flipping incoming 
or outgoing indices: 

(-Li)klmn = \^J K^iJklnm = \^) K^iJlkmn • (^6) 

As a consequence one finds also a simple transformation property of the re- 
duced amplitudes, a\ (u,t) = (— l) l+I a\ (t,u), under the interchange of in- 
coming or the outgoing particles. Furthermore the tensors Tj are chosen such 
that the Tj>3 are kinematically suppressed by one power of u or t if multiplied 
with on-shell spinors. In appendix A the invariant amplitudes a\ (t,u) are 
uniquely reconstructed in terms of the nucleon phase shifts. 

In this work we focus on the isospin one channel and evaluate the a\ (u,t)- 
amplitudes to leading orders in our chiral expansion scheme. Results on the 
isospin zero amplitudes will be reported on in a forthcoming publication. Here 
we elaborate in particular on the s-wave partial wave amplitude Tri,5 ](p) which 
is parameterized uniquely in terms of the s-wave phase shift 5ng ](p). The s- 
wave amplitude can be constructed in terms of our reduced amplitudes 



Tm<p)=l £ T mte t ) = $.Tj(?" fv -0 < 57 > 



where 
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/ \ (1) / \ ^ — t (1) I \ u ~^~ t (1) I \ I \ 

T[i So ] [u, t) = a\' (u, t) a\ (u, t) a\ (u, t) (58) 

and t = -2p 2 (l -cos 9), u = -2p 2 (l + cos9) (see (A.5,A.8)). 

We proceed and derive the leading terms of the scattering amplitude (39). The 
central interaction vertex relevant for our work is the pion-nucleon coupling 
vertex extracted from the chiral Lagrangian 



C = g A \N lblll (d^-r)N= 9 -^Ni 1 ^-TN + --- (59) 

Jit 

where the coupling constant A is identified with the pion decay constant A^ 1 = 
2 /„. ~ 186 MeV. For the axial vector coupling constant of the nucleon we take 
qa — 1.26. To the leading orders we will be working here the interaction kernel 
K n of (32) is identified with the one-pion exchange contribution 



K * = ~s~f2 ~~TZ — ~ V 1 - 3 ^oj + exchange term (60) 

° J it t — m^ 

following from (59) and Q = p\ — p[. The on-shell irreducible kernel K n of 
(44) follows with 



r> g 2 A ™ 2 75^75 (y %T> 

<J ] "' 2 {2 (ml - |) I 3fir (t) + I^{t) - 2 I 2 (m 2 )) (T\ + T 2 ) (p 1 + 9 V 1 



16 f, 
11 g\m 2 

384/4 



7i, w (Ti + T 2 ) (Vi + 9 Pi) + O (Q 2 ) + exchange terms (61) 



At leading orders it is given by the pion-nucleon vertex with pseudo-scalar 
coupling and contributions from scalar one-loop functions Iz^{t),l2,-rr{t),Ii,-n 
and ^(m 2 ) defined in appendix B. All we would like to point out here is that 
we find: 



(X + K n G K w ) | on _ shell = K v + K G K + O (Q 2 ) (62) 

by construction. Our result (61) clearly exemplifies the merit of transforming 
to the on-shell irreducible kernel: it permits to move the divergent pionic 
tadpole contribution J l7r into the local interaction vertex K a of (32) at the level 
of the interaction kernel. An interesting observation worked out in detail in 
appendix B.l concerns the chiral power of the pionic tadpole diagram 7 1;7r . As 
is demonstrated explicitly a consistent 1/m expansion implies ii ;7r ~ Q rather 
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than I\^ ~ Q 2 as one might expect naively from (3). Here this does not cause 
any complication since I\^ is absorbed into the bare coupling function. Note 
a further subtle property of the on-shell irreducible kernel K n : the scalar pion 
bubble loop l2,n{t) in (61), carries chiral power zero (see (B.10)). This does not 
contradict the chiral expansion since this term is canceled by a corresponding 
term in K n G K n by construction. 

The kernel K a was already given in (15) in terms of the energy dependent 
coupling function g(s). It is on-shell irreducible by construction. We further 
need to specify the subtraction interaction K$ of (36) and Kg of (41) 

K s=yt =c(s)g 2 A \ 2 T 1 V 1 , K s = c(s)g 2 A \ 2 T 1 V 1 , 

K a , s = g*,s(s) V^ = (<?(*) " c ( s ) 9a A 2 ) T i Pi 



(63) 



appropriate for s-wave scattering in the spin singlet channel. Here we allow for 
a s-dependence of the subtraction function c(s). The subtraction function c(s) 
is on the one hand required to justify our previously introduced L-counting 
rule and on the other hand ultimately leads to unitary results for the scattering 
amplitude. 

The leading orders as they contribute to our amplitudes are readily derived. 
The amplitude 7^ g(s) of (38) receives at chiral order Q contributions 



T 




TT,S- 



+ exchange terms (64) 

from the subtraction vertex Kg, the one-pion exchange, the pion nucleon box 
and triangle diagrams induced by Kg- In (64) the wavy line represents the pion 
exchange K w of (60) rather than K n of (61) for convenience. Transforming back 
to the on-shell irreducible kernel K n is straightforward and only required for 
the renormalization process. The full scattering amplitude follows now from 
(39). We note that all reduced amplitudes a\ (u,t) receive contributions from 
T n> s whereas the amplitude T a ^ affects only Oi (u,t): 

a?\u,t) = 2g 2 A \ 2 c(s)(l + 2g 2 A \ 2 VAs))+2g A \ i c 2 (s)J(s) 
+ g 4 A \ 4 (B 1 (u,t) + B 1 (t,u) 
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+ (l + g 2 A X 2 V^si-s) + g A A 4 V^ s (s)) T^s) , 
a§-\u,t)=g\\* (B 2 (u,t) - B 2 (t,u)) , 



.W 



2 2 \2 



a 3 '(u,t)=m g A X 



1 



t — ml u — ml 



+ g A \UB 3 (u,t) + B 3 (t,u)), 



,(!) 



2 „2 \2 



o4 (u,t) =m g A X 



1 



1 



, , j+^A 4 5 4 Kt)-5 4 (t )U ) , 

u — mi t — mi x ' 



a ( i\u,t)=m 2 g A X 



2 2 \2 



1 



t — mi 



+ 



u — mi 



+ g A X 4 (B 5 (u,t) + B 5 (t 



u , 



(65) 



The functions Bi(u,t) represent the contribution from the pion-nucleon box 
diagram in (64). The precise definition for all loop functions can be found in 
appendix B. The one-pion exchange contribution is made explicit in (65) and 
confirms the transformation property of the amplitudes a,i(t,u) under inter- 
change of t and u. The vertex loop functions Kr",s( s ) i n (65) are introduced 
with respect to the interaction kernel K^ t s of (36). The one and two- loop 
vertex functions K-,s(s) and K-^s^s) are then given by 



K,s(s) 



,2 \2 



gU 

V n (s) + c(s) J(s) 



X'XK 



V^As) 



<?1A 4 



+ 



X 



+ 





V^(s) + c(s) (J(s) V w (s) + .Us)) + c 2 (s) J 2 (s) 



(66) 



in terms of the vertex functions V w (s) and V n7T (s) defined with respect to the 
unsubtracted interaction kernel K n of (60), the not yet specified subtraction 
function c(s) and the one-loop bubble function J(s) introduced in (16). The 
pole term in (65) 



T*M 



(67) 



9a,s( s "> ~ J ( s ) ~ 9 a A2 J ir,s(s) - g\ A 4 J^, s {s) + ... 
is determined by the bubble functions J n ™,s(s) introduced with respect to 
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the interaction kernel K n ,s of (36). Again the subtracted two and three- loop 
bubble integrals Jn,s(s) and Jnir,s( s ) can be expressed as follows 



J«A s ) 



j w (s) + c(s)J' z (s 



)o<*)oo< 



J WW .,<?(S) 



g\^ \ 



+ 



)oo( * )ooo( 



= ^(s) + 2 c(s) J{s) J v (s) + c 2 (s) J 3 (s) (68) 

where J n (s) and J ffff (s) are the 'unsubtracted' loop functions. 

We will demonstrate that our L-counting rule can be applied for the loop 
functions 



V^M ~ Q n > >W S ) - -W^o) - Q" +1 (69) 

provided that c = c with c = 1/2 + (Q 2 ) and the subtraction point so ~ 4m 2 
sufficiently close to threshold. According to the L-counting rule the leading 
chiral power of the vertex functions K-™,s(s) and bubble functions J n ™,s( s ) — 
Jn™,s(so) is given by the number of independent loops. We return to the chiral 
powers of these loop functions in more detail below when we present explicit 
analytic results for all loop functions introduced here. For more details on the 
evaluation of the loop functions we refer the reader to appendix B. 

We proceed and systematically rewrite the bare loop functions J(p), Jn,s(p) 
and V^sip) of (65) and (67) in terms of renormalized loop functions Jr(p, z ), 
Jit,r(P' z o) an d K-,r(p)- Note that from now on we consider all loop functions 
as a function of the small momentum p rather than the Mandelstam variable 
s = 4 (m 2 +p 2 ). For simplicity we consider in this section only the leading term 
in the 1/m expansion and systematically drop correction terms. Technical 
details on the regularization scheme and on sub leading terms in the 1/m 
expansion can be found in appendix B. We find: 

J(p) = j(-i Zo ) + j R fa Zo ) , V VtS {p) = V VtR (p) + (c - |) J(p) , 
J*,s(p) = j *APi *o) + M-izo) + (c - ±) J 2 (p) . (70) 
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Note that here we also drop a contribution from the pion production cut in 
J-n,s{p) ( see (B.40)). The subtraction point z will ultimately be identified with 
the pseudo-bound state pole position in the chiral limit. If z > the subtrac- 
tion point lies on the second Riemann sheet of p = v/s/4 — m 2 . The renor- 
malized bubble loop functions Jr{j>, Zq) and Jtt,r(p, zq) vanish at p = — i zq by 
definition. The analytic expressions for the renormalized loop functions read: 



Jr{P,zo) = -—[ - zo + ip) , V^ R (p) = — 1 In [l-i — 

J*,r(P, zo) = -2 . r 2 In — - . 71 

All divergencies sit in the one-loop and two-loop subtraction coefficients J(—i Zq) 
and J n (—i z ) which may be estimated in terms of a cutoff parameter A intro- 
duced in the s-channel dispersion relation (see (17)): 



J(-i zo) - - A , U-i zo) - ^$ In (^] . (72) 

We point out that the vertex function V^^sip) is rendered finite by the unique 
choice c = 1/2 whereas the two- loop bubble J Wt s(p) requires a further sub- 
traction even at c = 1/2. This reflects our renormalization procedure (41) and 
(42). The vertex function is renormalized by an appropriate 'counter loop' 
leading to the renormalization point independent result K-,.r(p)- The bub- 
ble function Jk,s(p) is renormalized in two steps. The counter loops lead to 
a function J n (p) which defines the renormalized loop modulo a subtraction 
Jtt,r(p, z o) — Jnip) ~ Jn(—i z o)- This leads to a subtraction point dependent 
result for J n ,R(p,zo)- Note that here we are implicitly assuming that the on- 
shell irreducible interaction kernel K n is applied (see appendix B.3). 

Consistency of our L-counting rule requires the subtraction constant c to be 
rather close to one half. For example if c deviates largely from 1/2 the term 
proportional to (c — |) J(— iz ) Jr(p,z ) in our expression for J WtS must not 
be counted as chiral order two since g\ A 2 J(— i z ) > 1. We count (c — c) ~ Q 
and therewith confirm Kr,s(p) ~ Q an d J n ,s(p) ~ J-n,s(—izo) ~ Q 2 ■ Note that 
according to (70) J n ,s{p) carries in fact the leading chiral power Q. This result 
generalizes to the n-th order bubble J n ™,s(p)- The anomalous contribution to 
all orders is readily anticipated 



J^,s(p) = {c- l 2 ) n J n+ \-^o) + 0(Q n+1 ) . (73) 

and leads to J^^p) ~ Q n - 
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We proceed and express the not yet specified two-loop function V nn< s(p) an d 
three-loop function J^^sip) °f (65) in terms of the renormalized loops K-7r,i?(p) 
and J n ^ R (p, z ) 



V^sip) = V^ R (p) + (c - i) J(p) V n , R (p) + (c - |) 2 J 2 (p) 

+ ( C ~~ 2) (-^r(-« 2o) + ^,H(P, ^0)) 
Jn7T,s(p) = JttttAPi Z o) + -A«r(-* «0) + (c - I) J 3 (P) 

+ 2 (c - 1) J(p) ( J^(-i 2q) + ^r,fl(p, *o)) (74) 

at leading order in the 1/ra expansion. We suppress here sub leading contribu- 
tions from pion production cuts. All divergencies in (74) sit in the previously 
introduced power divergent one-loop bubble subtraction J(— i z ) and the log- 
divergent two-loop bubble subtraction J 7r (— iz ). The three-loop subtraction 
constant J wn (—iz ) is in fact finite: 



J ""<-^f#('"< 16 > + ©)- < 75 > 

The multiple powers of J n (—iz ) ~ A n in (74) reflect the presence of over- 
lapping divergencies in the two-loop vertex and three-loop bubble functions. 
They are easily disentangled by proper loop subtractions and as shown in ap- 
pendix B. Note that our result (74) confirms the anomalous part of J^^sip) 
as anticipated in (73). It is straightforward to derive analytic expressions for 
the renormalized loop functions. Again we find it most economic to deduce 
the analytic form of the loop functions from their imaginary parts by means 
of their known analytic structure. We drop 1/m correction terms and present 
the leading renormalized loop functions 



. . m 2 m? m? ( / 1 ip\ / 1 m* 



(47r) 2 4p 2 \ \2 m 7r J \2m 7r — ip 



1 \ 12/ m -K — '2ip\ , ( 1 ip 



+*te^ K)K<»> (76) 
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with J-Kn^ip, z) = Jmrip) — Jmr(~i z) 7 . Again we confirm the anticipated chi- 
ral powers of the loop functions with V n7T) s(p) ~ Q 2 and J mX! s{p)~Jir-n,s{~i z o) rv - 
Q 3 . The analytic results (76) explicitly demonstrate the presence of a branch 
point at p — —im^ in the two-loop vertex function V n7! -,ii(p) and three-loop 
bubble function Jmr,R(p,zo)- Note that the polylog function Li 2 (x) has one 
branch point at x — 1. It remains the contribution from the pion-nucleon box 
diagram to the s-wave scattering amplitude. The relevant combination 



u — t u + t _ 1 

— — B A (u, t) - -— B 5 (u, t) = B(u, t) - V n , R (p) + - 
Am z Am z 4 



is conveniently expressed in terms of the reduced pion-nucleon box function 
B(u, t), the renormalized vertex loop V n ^(p) and J(p). We derive the analytic 
expression 



t-./ \ 1 mrnt ( . I t(u + t) 

B(u,t) — - ; =\i arctanh \, —- — 

1 } 47r^_ t6(M )V \\\ Ku,t) 



' y/-tb(u,t) __ y/-u-t ( 



1 
H — arctan 

2 \ Ami m„ 




+ - arctan ( JE^M + ^^ ( x _ 

-f- — dlCldll I - - | 1 — 

2 I 4 m; m w \ 

b{u,t) = 4ml-4ml{u + t)+t{u + t) . (78) 

The result (78) confirms the expected chiral order of the reduced diagram with 
B(u,t)~Q. 

Finally we collect all terms introduced in (65) and (67) relevant for s-wave 
scattering according to (58) 



T [1So] (u,t)=2g 2 A X 2 (c- \) (l + (c-\) J(p))+4g A \ 4 (c - \) V n , R (p) 

-\s\\ 2 ( -^ + r^M + A A 4 (B(u, t) + B(t, u 
2 \u — m z t — ml) v 

(l + g\ A 2 V n , s (p)) 2 + 9a A 2 V mtS (p) 

4- 9 - - ( 7 C )) 

g;Us)-J(p)-g 2 A X 2 J^ s (p)-9 4 A^J^s(p) 



7 Here we recall the definition of the polylog function 

T . , , } , ln(l -t) ^x k 
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at chiral order Q. According to the discussion in section 2 we refrain from 
expanding the denominator in (79). Expression (79) serves as a convenient 
starting point for any further expansion implied by scenario I and II of section 
2. 

It is worth pointing out a generic property of the renormalized loop functions 
Vtt",r(p) and Jn",R.(p,Zo). Since they satisfy a dispersion relation in the Man- 
delstam variable s they are in fact analytic in p except for multiple branch 
points at p = — i km n /2 with > n > k . This convenient property can be ap- 
plied to derive the analytic form of the loop function in terms of its imaginary 
part (see appendix B). The s-wave scattering amplitude Tns ](p) shows also 
branch points on the upper complex half plane with p = +ikm n /2. These 
branch points are, however, exclusively induced by the s-wave projected mul- 
tiple pion exchanges of T nj s- We perform the angle average of the one-pion 
exchange contribution and B(u,t) according to (57) 



*W = -7/^lT^ + 7^l=g>"(l + g 



r ail 

Rnn{p)=J J^B{ U ,t) 




mm^ m 3 






2 to^ + i p 

+21n f^L±il!lU arcta n 

\ 4 m^ + Ap z J \ \ m Tl 

+2 i ln(2) In ^ + ^ + 2 i ln 2 (2) + i^j . (80) 

The expressions (80) confirm that indeed the function R n (p) shows branch 
points at p — ±im 7T /2 and R nn (p) branch points at p — ±im 7T and p = 
±i 171^/2 as expected. 

We conclude that the chiral L-counting rule derived by means of the s-channel 
dispersion relation in section 2 is confirmed explicitly for the loop functions 
considered here. Let us shed some light on the convergence property of the 
chiral expansion. A convenient measure is the ratio of the one-to-two loop and 
two-to-three loop subtraction constants: 



2 ,2 ^,s(0) 2 2 mm n m w [A 2 

9a A ' ' ^ a A A — — In — 

A J(0) A 4vr A U 
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2 ,2 J**,S( Q ) 2 ,2 mm n In (4) 






at c = c. The typical cutoff A ~ 600 MeV implies reasonable small ratios of 
about 0.27 and 0.35 for the leading and subleading ratio respectively. Finally 
consider also the ratio of leading to subleading vertex functions: 



g2 A 2 ^§ = g\*^ In (2) c 0.28 (82) 

evaluated at p 2 = 0. All together we find good convergence behavior of the 
chiral expansion. The underlying expansion parameter appears to be roughly 
1/3. Before presenting further details of our chiral expansion scheme we com- 
pare the analytic results with a calculation by Fogel[29] obtained in the static 
pion limit. 



3. 1 Static pion limit 



It is instructive to compare our result (65) with a corresponding calculation 
performed with static pion propagators. As an example we consider here the 
one-pion exchange contributions R w (p) and R wn (p) (see (80)). The limit of 
static pions is accessed most economically by considering the Schrodinger 
equation of a Yukawa potential 



il>"(r) + (p 2 + g n m n ——) ^(r) = . (83) 

with the coupling constant g^ = g\m 7T m/(167r f 2 ) ~ 0.45. This problem was 
solved analytically to lowest orders in the coupling constant by Fogel [29] . The 
s-wave scattering phase shift 5[is ]{p) 



»/(p) 



tan5 [l5o] (p) = ^4^( (84) 



is expressed in terms of the Jost function, f(p), which is the solution of (83) 
with asymptotic behavior ~ e~ ipr evaluated at r = 0. In [29] the Jost function 
is evaluated perturbatively 



f(p) = i + J2g:fn( P ) (85) 

n=\ 

with analytic results given for leading and subleading terms: 
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hip =^ln 1 + 2— , 
2 p V m^J 



2p 

, i / 1 / ,• , , x tj-2 (2ip — m n 



/ 2 (p) = -T^U ln U + 2— +T7T- Li 




±2 



12 " \2ip + m^ J 

(86) 

It is found [29,30] that for physical coupling g„- ~ 0.45 the perturbative expan- 
sion of the Jost function is rapidly converging for momenta 2 p > g n m n ~ 70 
MeV. The convergence property is in fact improving significantly as p gets 
larger. The next moment fzip) was evaluated by Friman [30]. 

We confront our result with the analytic expressions of Fogel by matching the 
s-wave scattering amplitude, Tng ](p), at leading and subleading order in the 
coupling constant g n : 



T ^ p) = V^Mf(p)-i*f(p) (87) 

We identify 



^(p)~^3/i(p), 
2p 

RUP) - ^ f9f/ 2 (p) " 9/i(p) (»/i(p) " i*fi(p)j\ (88) 

and conclude that, as to be expected, the terms agree identically. For a discus- 
sion and evaluation of 1/m correction terms as they evolve from relativistic 
Feynman diagrams see appendix B.l. 



3.2 Expansion of effective coupling function 



In this section we will further prepare the ground to generalize (25) for the 
case in which the pion dynamics is included explicitly in the theory. Starting 
from (79) we will specify how to identify and expand the renormalized coupling 
function. 

The scattering amplitude (79) involves the subtracted coupling function g a s = 
g — g\\ 2 c. Within our renormalization procedure (41) and (42), which holds 
for both scenarios I) and II), the loop subtraction coefficient c is determined 
by enforcing finite vertex functions Kr«,s(p)- At leading orders this is achieved 
with c = 1/2 + 0(Q 2 ) corresponding to the 'local' part of the one-pion ex- 
change potential (see (79)). The remaining divergencies sit in the bubble loop 
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functions J n n ,s{s) (see (70,74)) and can be absorbed into the bare coupling 
function g{p^m 7T ) by a subtraction procedure. To this end it is convenient to 
introduce the auxiliary coupling function g H (s, m n ] c), the basic building block 
in (79), 

5- H 1 (s,m 7r ;c) = g'sfam*; c) - J(-iz ) 

-g\ A 2 J Vj s(-i 2b) - 9a ^ 4 J™,s(-i z ) H , 

9r(p, rn n ) = 9h (s, m^ c) (89) 

which then defines the renormalized coupling function if evaluated at c = c. 
In (89) we consider the bubble loops J^ s {ji) as a function of the momentum 
p. The auxiliary coupling function g H (s,m n ;c) depends on the subtraction c 
via the coupling function g ff s and via the explicit dependence of the bubble 
loop functions J n n,s(p)- The subtraction point z = z(m n = 0) is identified 
with the pseudo-bound state pole position in the chiral limit. 

The bare coupling function g(s, m 2 ) in gn{s, m n ; c) permits a Taylor expansion 
in p 2 and m 2 . In scenario I) the bare coupling function g(s,m%) and the 
renormalized coupling function gR(p,m n ) are expanded as follows: 

2\_ \2 V^ fl (2fc) / 2 | Jl\ n ^2k 



g(s,ml) = \*j:ttn } (p 2 + z l 



m„ 



"0 1 ""-re > 
n,fc 

9£<P,m«) = £ E CS (P 2 + <f m* . (90) 

n,k 

Note that in (90) we slightly reorganized the expansion of the bare cou- 
pling function as compared to (15). The function g{s^m 1T ) represents an in- 
finite hierarchy of interaction terms of the chiral Lagrangian with increas- 
ing number of derivative couplings and increasing degree of explicit chiral 
symmetry breaking. The naturalness assumption leads to fiffi ~ A^at" 1 an d 



(-0 



A na t. = A„ at ~ rrip/2 at given cutoff A ~ A 



rial. 



The renormalized coupling constants Q 2 n can be expressed in terms of the bare 
parameters (3 2n an d the subtraction coefficients J^n : s{— izq). We expand the 
inverse auxiliary coupling function gj 1 1 (s,m n ;c) around the point c = c (see 
(48)). The leading terms read: 



ti 2) < + $> (p 2 + zl) 



g(s, m n ) -cg\X> $) _ g y 2 (/#> - g y 2 ) 



+ 



9a(c-c) 9a\c- c 



2 

2 
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2 Ale /"i^^^) |C W (91, 



where we count c — c ~ Q. The divergent subtraction coefficients can now be 
absorbed into the bare parameters fi\ n . The linear divergence in J(—izo) is 
absorbed into the bare coupling fy and the log-divergence in J 7r (— i z ) is ab- 
sorbed into fy ( see (72)). The leading order renormalized coupling constants 
Co i C2 an d Co can now be identified 



A 2 J(-iz ) = T - M — ^ , C: 



(0) _ 2n(3J 0) /X 2 

2 / _mi „ ,_\2 ' 



(/^-si/2) ' " -(4 0) -.i/2) 



d 2) = -^ A 2 =- 



2n U-iz ) 2tt/35 2) /A 2 



m mH. /~rm o ,„\ 2 



m (4°' - Aft 



tf^-g^HMpl, (92 ) 

7T 

Note that the renormalized coupling function gjt{p, m n ) includes non-analytic 
terms in the squared pion mass (e.g. m^Q ) from the bubble loop sub- 
traction coefficients J n ",s(— izo) not present in the bare coupling function 
g(s,m 2 ). Moreover we point out that the inverse auxiliary coupling function 
g]j (s,m 7T ; c) and the inverse renormalized coupling function g^ (p,m n ) agree 
up to chiral order Q 2 : 



4» ™ (J> , ,A j- ™ A 2 ) ^2 , n f n s 



9h\s, m^ c) = (?> — (p 2 + 4) + — Co' K + O (Q* . (93) 

Z7T V '111 v/ 

This follows since (92) implies that the power divergent terms proportional 
to (c — c) n with n — 1, 2 in (91) are canceled by matching contributions from 
J-K,s(—i z o) an d J-K-K,s{~i z o) ( see (70,73)). For the generalization of (93), which 
does show an explicit dependence on c—c, valid up to chiral order Q 3 see (C.5). 

We turn to scenario II). Most economically it may be considered as a particular 
resummation of scenario I) with: 



CO oo 

£ C£ < = -£, {£ + zl) ~ n + £ #2 mj (94) 

fc=0 fc=0 
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where \i a = Hcr(m n ) and £ CT = $, ff {m^). The bare coupling function g(s,m 2 ) 
and the renormalized inverse coupling functions gf j l II {p^m ll ) of scenario II) 
are expanded around p 2 = fi 2 . At given order (p 2 — fi 2 ) n the terms can then 
always be regrouped according to (94). Note that the form of the renormal- 
ized coupling functions in (90,94) avoids a fine tuning of parameters which in 
(25) is required for the formation of the pseudo-bound state pole. The param- 
eter \x a of scenario II) carries physical significance since it will describe the 
empirical fact that s-wave scattering at relative momentum p ~ 350 MeV is 
interaction free. This phenomena reflects the well established subtle cancella- 
tion of attractive cr-exchange and repulsive u;-exchange in phenomenological 
nucleon-nucleon potentials [24,25]. The expression (94) clearly demonstrate 
that as the value of /i CT decreases the effective chiral power of the first term in 
gf^n changes dramatically: for /i CT ^> A nat we find Q 2 whereas for fi a <C A nat 
we are led to Q°. Note however that at momenta p 2 <C fi 2 the 'leading' term is 
suppressed by the factor z 2 /fj, 2 and therefore the 'subleading' term, in particu- 

(2) 

lar the Q -term, becomes equally important. Recall that we expect /i CT ~ 350 
MeV and zq ~ z ~ 8 MeV. Thus an expansion in ratios of 'subleading' over 
'leading' term is converging for momenta p 2 ~ fi 2 but turns useless for mo- 
menta p 2 <C /i^. Our expansion for g^j, in contrast, is expected to converge 

well for momenta /i^ — A; ia l < p 2 < n 2 + A^l with A^J ~ 600 MeV (see sec- 
tion 2.2). This is further support for our claim that T^s must not be expanded 
in powers of T^s ( see (38)). 



3.3 Inefficient chiral expansion (CE) scheme 



In this section we develop the scenario I) in the presence of perturbative 
pion dynamics. The denominator of (67) is systematically expanded in chiral 
powers. We call such an approach 'CE' scheme for future reference. The mo- 
tivation of this section is twofold. On the one hand we show the failure of the 
CE scheme on the other hand we present our renormalization scheme at work. 
We will systematically express the scattering amplitude at subleading orders 
in terms of the physical parameters z , Q introduced in the previous section 
(see (97,92)) rather than in terms of the cutoff and bare coupling parameters. 
Therewith manifest scheme independent results are obtained. In particular we 
will demonstrate the independence of the scattering amplitude on the sub- 
traction function c(s) in a perturbative sense: if the scattering amplitude is 
expanded up to chiral order Q n with c — c ~ Q any residual c-dependence 
contributes at order Q n+l . The perturbative cancellation of any c-dependence 
demonstrates that physical results do not depend on how to distribute strength 
among K n> s and K a s in (36,37). For the full solution of the Bethe-Salpeter 
equation this is of course true by construction. For the truncated scattering 
amplitude this is, however, a non-trivial constraint. 
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According to (93) the auxiliary coupling function gJ £ 1 (s,m 7T ] c) carries leading 
chiral order two. This leads to an s-wave scattering amplitude of leading chiral 
order minus one: 



T m {p) = "TT^ + O (Q ) = — -V- + O (Q ) . (95) 

Jr{p,z ) v ' m p + izq V / 

The result (95) is certainly not a surprise. One recovers the effective range 
theory with r = 0. The pole in (95) at p — —izo is the regularization scheme 
dependent result of fine-tuned bare coupling strength and intrinsic cutoff pa- 
rameter. So far the chiral dynamics shows no quantitative predictive power. 
First non-trivial consequence are expected at subleading order. Note that the 
leading order result does not yet reveal the necessity to expand around c = c 
with c— 1/2 + 0(Q 2 ). This can first be seen at subleading order. 

We proceed and turn to the zeroth order amplitude. At this order the auxiliary 
coupling function gn(s, fn n ; c) can be identified with the renormalized coupling 
function gR{p,m n ) (see (93)). We collect all terms of chiral order zero: 



* A *{ l -<W5 + **&) +ov »- (96) 

The result (96) is expressed in terms of the physical parameters in gR{p,m n ) 
and a set of renormalized loop functions given in (71,80). We point out that 
in (96) the explicit c-dependence induced by the loop functions V^^ip) an d 
Jir,s(p) (see (70)) cancels identically upon expansion. Therefore the zeroth 
order scattering amplitude (96) is manifestly independent of the subtraction 
procedure. From the analysis of the effective range theory we expect the scat- 
tering amplitude to exhibit a single pole at p = — i z rather than at p = —iz . 
We conjecture an expansion of the pseudo-bound state pole position z(m n ) in 
powers of the pion mass: 



z(m 7T ) = z + ^ ( m m™ (97) 

m = 2 

with C, m ~ A~^ +1 and z ~ m n ~ Q. Therefore an expansion scheme which 
preserves the proper low energy structure of the scattering amplitude (27) 
requires a particular summation. In order to comply with (27) we expand all 
bare poles in terms of powers of the physical pole term and then impose the 
cancellation of all multiple poles at p = — i z. The expansion of the leading 
order pole reads: 
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J R (p,z ) Jr{p,z) ~ 2 Ti J R (p,z) " '' 2n 2 J R (p,z 

m m^. 
2 it J R (p,z 



+ 2 Cs — 127^ + --- (^ 



where we apply our notation (97). At subleading order this leads to ( 2 = Co • 

We now turn to the comparison of our result (96) with the effective range 
theory. It is instructive to make the analytic structure of our result more 
explicit. We rewrite (96) and present the central result of this section in terms 
of the reduced functions R V7r and R Jw and the one-pion exchange contribution 
R^p) (see (80)): 



T[i So ] (p) = — — + — Q 

m p + iz m 

+ 2g\\ 2 (R^p) - 2 RW(p, zo) + R%{p, z )) + O (Q) . (99) 

Here we introduced a convenient notation: for any given loop function F(p, z) 
we define the n-th order residual function, Rp (p,z), by: 



n— 1 1 

RP (p, z) = J R n (p, z) (F(p, z )-J2 V ,(p + iz) k d k p F(-i z, z)) (100) 

fc=0 

Note that for a given loop function F(p, z) which is regular at p — —iz also 
the set of residual functions Rp (p, z) is regular at p = — i z by construction. 
The pole residuum w is: 



22/ / \ 27r d . s 

w = l + g A \ [2V 7r>R {-izo) -z-J v ,R{.-izo,zo, 

v m op 

-2(i 0) z + O(Q 2 ) . (101) 

With (99) we arrive at a representation of the scattering amplitude anticipated 
in the previous section when discussing the low energy limit. The amplitude 
is the sum of the pseudo-bound state pole and a smooth non-polynomial re- 
mainder, which takes care of the non-trivial pion dynamics. 

In Fig. 1) we show the result of a low energy fit to the empirical scattering 
amplitude. For momenta p < 50 MeV the scattering amplitude is well re- 
produced. It is also shown that the unitarity condition is violated strongly 
at momenta p < 150 MeV. This is a disappointing result since already the 
effective range approximation works quite accurately up to momenta p ~ 150 
MeV. We conclude that it would be misleading to confront the CE result (99) 
with the empirical scattering phase shift. 
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Fig. 1. The left hand figure shows the real part of the np-scattering amplitude 
f[iS ](j>) = VsT[i So }(p)/(8tt). The right hand figure shows the imaginary part of 
the np-scattering amplitude Q/ri5 i(p) and the unitarity function —^sfn s Ap)/p, 
which is one if unitarity holds. The solid lines are from the PWA93 solution of [4] 
and the dashed lines are the result of the zeroth order CE scheme with Q 
fm, z = 8.35 MeV and z = 7.89 MeV in (99). 



0.380 



It is instructive to analyze our result in more detail and apply a further nu- 
merical expansion to the pole residuum w and the effective scattering length 
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+ 0(Q) (102) 



The free parameter Q and zq are now adjusted to reproduce the empirical 
residuum w ~ 0.90 and the effective scattering length b e g ~ —1.15 fm (see 
(29)). We obtain the values Q — 0-38 fm and z ~ 8.35 MeV. The scheme 
dependent bare parameters /3q , fy ano ^ Po follow with (92). The remaining 



.(2) 



(2) 



physical parameter Q follows with z = z + Q m i + 

Q ~ —0.0047 fm. It is important to inspect also the next moment 



which leads to 
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Fig. 2. The left hand figure shows real part of the np-scattering amplitude 
f[iS ](j>) = VsT[i So }(p)/(8tt). The right hand figure shows imaginary part of the 
np-scattering amplitude Q/rig i(p) and the unitarity function —^sfn s \{p)/p, which 
is one if unitarity holds. The solid lines are from the PWA93 solution of [4] and 
the dashed lines are the result of the first order CE scheme with Q 
^ 2 ) = 1.771 f m and z = -45.68 MeV (see (C.8)). 



-1.363 fm, 



b (1) 
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(103) 



and compare it with the result from effective range theory (see (29)). We 
conclude that (103) predicts the wrong sign of this term and therewith confirm 
our result that (99) is in strong conflict with the unitarity condition 



^ J I 1 So] 



(P) 



■-$J R (p,Zo) 



(104) 



In appendix C the reader may find details of the NNLO terms in the CE 
scheme. In Fig 2 we present the result of a low energy fit. The description of 
the scattering amplitude is somewhat improved for p < 100 MeV as compared 
to the zeroth order result. In particular the unitarity condition is met to better 
accuracy. However the first order calculation clearly signals the breakdown of 
the CE scheme at momenta p ~ 150 MeV. To be explicit: at the chiral orders 
we are working here (see (65)) the strength of the effective pole at p = iz 
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(see (30)) must be represented by the pion-nucleon box function R nn (p) (see 
(80)) since it is the only function showing a branch point in the upper complex 
half plane at p — % m n ~ iz. Recall here that all vertex functions and bubble 
functions have branch cuts only on the lower complex half plane. Since the 
strength of the cut in R n7T (p) is much too weak to represent the strength of the 
effective pole at p = iz the CE scheme fails for p<z~ 150 MeV. Note also 
that the parameters Q an d zq of the zeroth order and first order calculation 
are strongly renormalized. This points at the poor convergence properties of 
the CE scheme even at low momenta p < 150 MeV. This could have been 
anticipated by inspecting the ratio of the renormalized functions 



g2 A2 <^M = A2 mp, + Q ^ Q g 
Jr(0,z) 2tt 

at p — 0. From our calculation we conclude that in fact we have to change 
scheme and reformulate our effective field theory along the lines already an- 
ticipated. We emphasize that the result (105) does not imply that successive 
unitary iterations of K Wy s need to be summed to all orders. As shown before 
here typical ratios are well behaved (see (81,82)). The ratio (105) is to be 
taken as further evidence that the properly renormalized denominator (67) 
must not be expanded. 

We close this section by a comment on the so called KSW scheme proposed 
by Kaplan, Savage and Wise [15]. In the KSW scheme the scattering phase 
Sns ](p) is expanded perturbatively in the pion dynamics imposing the correct 
pseudo-bound state structure. It followed from our CE scheme upon a further 
expansion. For independent calculations in the KSW scheme at subleading 
orders see [31]. In our work we insist on the CE scheme since the unitar- 
ity condition (104) provides a powerful consistency check for the CE scheme 
which is obscured in the KSW scheme. If the CE scheme showed good con- 
vergence properties it would be equivalent to the KSW scheme. The rather 
poor convergence properties of the CE scheme necessarily implies that also 
the KSW scheme is poorly convergent. We conclude therefore that neither the 
CE nor the KSW scheme provides an improvement upon the effective range 
theory. Furthermore note that by expanding the phase shift directly rather 
than the amplitude makes it more difficult to detect a potential breakdown of 
the scheme. This is clearly demonstrated by the effective range theory. Here 
the second low momentum pole at p = iz (see (27,30)) restricts the applica- 
bility of the CE and the KSW scheme to momenta p < \z\. The breakdown 
is triggered by a 'weak' cut singularity in the KSW scheme as compared to a 
'strong' pole singularity in the CE scheme: 




<Vs ](p) = 7T- arctan + o r H • ( 106 ) 
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In the KSW scheme the branch cut of the arctan function at ± % leads to the 
convergence bound p < \z\ — |2/r| + 0(1/ a). We conclude that in order to 
uncover the physical significance of non-perturbative strength at p ~ % z it is 
advantageous to choose the CE scheme. 



3.4 Unitary chiral expansion (UCE) scheme 



In the previous sections we collected mounting evidence that the denominator 
in (38) or (79) must not be expanded further. We have demonstrated that the 
CE scheme, based on such an expansion, is not applicable in the spin singlet 
channel due to poor convergence properties. In this central section of this work 
our new scheme already outlined in section 2.4 is constructed explicitly. It 
shows much improved convergence properties and comprises unitarity exactly. 
Formally the scheme arises if the renormalized coupling, gR 1 (p,m- K ), changes 
according to (94) its chiral power from plus two at very small momenta p <C 
p, a to zero at intermediate small momenta p ~ fi a . As a consequence the 
denominator in (79) must not be expanded if the effective field theory is to 
be applied for momenta p < fi a + A nat .. Keeping the denominator on the 
one hand takes care of the proper chiral power of the renormalized coupling 
at intermediate small momenta p ~ yU CT and on the other hand induces a 
particular resummation for the low momentum regime which greatly improves 
its convergence property. We refer to this approach as the 'UCE'-scheme since, 
as we will elaborate on in detail below, it in addition leads to exact unitarity. 

The UCE scheme should be viewed as a systematic summation of the CE 
scheme (see (94)). The renormalization procedure is identical to the renormal- 
ization procedure of the CE scheme: all vertex functions V 7T n tS (p) are renormal- 
ized by appropriate 'loop' subtractions induced by a properly constructed local 
subtraction kernel Kg. We emphasize again that this procedure leads to renor- 
malization scale independent results for the nth order vertex functions V n n R (p) 
and the nth order box functions R^iji). The remaining divergencies of the nth 
order bubble functions J 7T «,s(p) are renormalized by a subtraction procedure 
at p = — i Zq. Thus the renormalized bubble functions Jn n ,R.(p, ^o) show an ex- 
plicit dependence on the subtraction point Zq or more generally they depend 
on some renormalization scale. At strictly Ks = Kg we obtain all together 
renormalization scale independent results since the renormalization scale can 
be absorbed into the renormalized coupling function g R (p, m n ). As advocated 
in section 2.4 it is advantageous to exploit the freedom to not strictly require 
Ks = Ks- For the CE scheme we have established that the scattering ampli- 
tude evaluated to a given chiral order is independent on (Ks— Kg) n if we count 
Ks — Ks ~ Q- In the UCE scheme this does not hold anymore since it is an in- 
finite resummation of the CE scheme. The residual dependence on (Kg — Ks) 71 
is fixed as to implement unitarity exactly at given order in the expansion. In 
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order not to spoil the renormalization procedure the residual dependence on 
(Kg — K s ) n 'acts' exclusively on finite parts as defined in (47) and (48). This 
amounts to dropping all divergent terms proportional to (c — 1/2) J(—i zq) in 
the loop functions (70). More explicitly we arrive at: 



T [lSo] (p) = 2g 2 A \ 2 (R n (p) + c-l) 

(107) 



2, 

2 + Ag 2 A X 2 V 1t , s (p : z ) 



9r(p, m n) - Jr(p, zo) - g\ A 2 J v ,s(p, zo) 
with the renormalized loop functions 



V^slp, z ) = (c-\) Jr(p, zq) + V KjR (p) , 

J«,S(P, Zo) = (c - |) J R (p, Zq) + J n ,R{P' zo) ■ (108) 

The renormalized coupling g R (p, m w ) is given by (94) and the loop functions in 
(71,80). Note that according to our definition (48) and (91) the finite 'part' of 
the auxiliary coupling function gH(s,m n ; c(s)) of (89) must be identified with 
the renormalized coupling gRiPjin^). The function c(s) is now unambiguously 
determined by the unitarity constraint. According to (49) the function c(s), 
representing a separable interaction, is constructed such that the s-wave and 
on-shell projection of T n s vanishes at the given order 



(T^ s )=2g 2 A X 2 {R w (p) + c-±)=0 + O(Q) . (109) 

We find that the unique choice satisfying (109) 



c( s ) = ±-R K (p) + 0(Q) (110) 

renders the scattering amplitude (107) unitary. Note that the subtraction func- 
tion c(s) — 1/2 now carries chiral order zero. This does not contradict the 
L-counting rule because the subtraction c(s) is now only active on convergent 
parts of the loops. We emphasize that the construction (110) holds irrespective 
of the form of the renormalized coupling g R (p, m n ). 

The mechanism which leads to a unitary scattering amplitude is readily un- 
derstood. The vertex loop function V^#(p), the bubble function J^r^Zo), 
and the one-pion exchange contribution R n (p) are not independent entities. 
They are interrelated by causality. Cutkosky's cutting rules [23] lead to the 
identities 



R n (p) 3 J R (p, zq) = 3 V K>R (p) 
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2(9W Wifl (p)) (QU R (p,z j)=QlJ n , R (p,z ) (111) 

which are consistent with our analytic results (71,80) and directly confirm 
the unitarity condition (104). As a consequence the phase shift can now be 
identified unambiguously 



, x /n ^n g R l (p,m ir )-g A \ 2 ?R:J 1 r,s(p,Zo) ( s 

p oa**p^(p) = -^ __^______ (112) 

where we introduced 



9r(p, m*) = g R (p, m w ) - 9fc J R (p, z ) (113) 

for notational convenience. 

We point out that any further expansion of (112) is necessarily in conflict 
with unitarity (or even worse with the complicated low energy structure of 
the scattering amplitude) since it would invalidate the derivation of our result 
(112). Finally it is instructive to consider the limit of a vanishing renormal- 
ized coupling function: g R (p,m 7T ) = leads to a zero phase shift S[is ](p) = 0. 
This demonstrates that the renormalized coupling function g R does not only 
represent the physics of short range repulsion (u;— exchange) and intermediate 
range attraction (a— exchange) it also encompasses some residual pion dy- 
namics. The important achievement lies in the fact that g R does not pick up 
pion induced singularities. In other words g R does not exhibit branch points 
at p — ±ni m n /2 with n = 1,2, 3, ... and therefore g R is expected to respect a 
Taylor expansion with the typical scale set by a rather large A nat . ~ 600 MeV. 

3.5 NNLO terms of UCE scheme 



We turn to the next terms in our expansion. Again all divergent terms in the 
loop functions (70,74) are dropped according to the argument given in the 
previous section. We derive 



T[i 5o] (p) = 2 g\ A 2 (R„(p) + c - \) + 2 g\ A 4 (c - \) J R (p, z ) 
+ 4^A 4 (c-\)v^ R {p) + 2g\\ A R^(p) 



2(l + g 2 A \ 2 V nt s(P,zo)) +4g A \ A V wn , s (P,zo) , , 

H =■ = (H4) 

9 R \p) - Jr(p, zq) - g\ A 2 J^^sip, z Q ) - g\ A 4 J^,s(p, z ) 



with the renormalized loop functions 
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Kr7r,s(p, zo) = V^nip) + (c - \j J 2 R (p, zq) 

+ (c - |) ( Jfl(p, Zo) V^flfc) + J,r,fl(p, Zo)) , 

Jmr,s(Pi z o) = J**,r(i>, zo) + ( c - I) ^|(p, zo) 

+ 2(c-i) J R (p, zo) J^p, z ) . (115) 



Before constructing the subtraction function c(s) at subleading order we col- 
lect the relevant causality identities predicted by Cutkosky's [23] cutting rules: 

^ Rmtip) = Rlip) $ Jr(p, Zo) , 

%V 7Tn , R (p) = ($tV 7T , R (p)) (3K, fi (p)) + (toRmip)) (SJr(p,zo)) , 
c *J™,R(p,z ) = 2(9tV mr , R (p)) (%J R (p,zo)) + \v n>R (p)\ 2 9 J*(p, z ) .(116) 



<y 



Again the identities (116) are respected by our analytic results (76,80). The 
subleading term of the subtraction function c(s) can now be constructed by 
the requirement 



(T^ s )=2g 2 A \ 2 (i? 7r (p) + c -I)+2^A 4 (c-\) J R (p, z ) 



4 g\ A 4 (c - \) V w , R (p) + 2 g A A 4 itU(p) 

+ (Q 2 ) (117) 



that the s-wave and on-shell projection of T n> s vanishes to the given order 
(49). The solution of (117) is now readily found with 



c(s) = 1 -- R„{p) -g\\ 2 ^ R^ip) - g\ A 2 R 2 n (p) U J R (p, z ) 

+ 2^A 2 ^(p)3?V;^) + 0(g 2 ) . (118) 

Explicit calculation indeed confirms that with (118) the s-wave scattering 
amplitude (114) satisfies unitarity. The phase shift can be identified unam- 
biguously and we finally arrive at the central result of our work: 



■ jt /n 4vr 9 R (p) ~9a ^ 2 &Jic,s(p, z o) -9a ^ & Jm,s(p, zo) 
P cot d[i So ] (p) = —j= — z T2 • 

V s ft (l + #2 A 2 ^^ ^)) + 2 ^4 A 4 gj K ^ s(p? Zq) 

(119) 

Note that in (119) the numerator has to be evaluated at chiral order three and 
the denominator at chiral order two. In other words in J n ,s{p) the subtraction 
function c(s) has to be included to chiral order one as given by (118) whereas 
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in Jnn,s(p) the subtraction function c(s) is to be included to chiral order zero 
as given by (110). 



4 Results and discussion 



We confront our unitary chiral expansion (UCE) with the partial wave anal- 
ysis (PWA93) of the Nijmegen group [4]. Our result at order Q° is given in 
(108,110,112) and at order Q 1 in (115,118,119). At leading chiral orders we 
encounter the renormalized coupling function 

characterized by five free parameters ^ CT ,Q >£o il^a and z . Here we absorb 
the finite contribution m% £q i n t° m "i £o f° r convenience. Our parameters 
are adjusted to reproduce the 1 5'o-partial wave in the pn-channel. With g n = 
g\mm- K l{V6'KJl) = 0.451, m n = 139.0 MeV and m = 938.9 MeV we perform 
a least square fit to the single energy 1 Sq partial wave as given in [4] with 
finite error bars. Tab. 1 presents the resulting parameters together with the 
X 2 -value of the fit. We observe that the x 2 °f the UCE result at order Q 1 is 
significantly improved upon the x 2 °f the UCE result at order Q° although 
there is no additional free parameter involved. The parameters are only weakly 
renormalized indicating good convergence properties and consistency with the 
naturalness assumption (see (26)). Moreover we confirm that the inclusion of 
pionic effects leads to a much better description of the partial wave. The x 2 is 
reduced by a factor of 5 from 25 to 5. The collum in Tab. 1 labeled 'no pions' 
refers to our calculation with g^ = 0. 

We emphasize that the parameter z has a physical interpretation. It represents 
the pseudo-bound state pole position in the chiral limit with m n = 0. The 
chiral limit leads in fact to a real bound state at p ~ i 37 MeV with binding 
energy e ~ Zq/itl ~ 1.5 MeV rather than a pseudo-bound state obtained with 
m n > 0. Our result demonstrates that the effective range parameter, r(}So), 
defined in (28), is driven by the physical scale p, a rather than the pion mass 
as one might expect naively. In the language of the a—,u— meson exchange 
picture we would say that the effective range is mainly determined by the 
subtle interplay of intermediate attraction and short range repulsion rather 
than by the pion exchange. The parameter \x a determines the zero of the 
phase shift and ( a 



5 m (p)^^-!^ (121) 
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the behavior close to p ~ \x a . In (121) we omitted pion induced terms for 
simplicity. We point out that the dynamic chiral power of our renormalized 
coupling function g R in (120) with g^ 1 ~ Q 2 for p <C ji a but g^ 1 ~ Q~ 2 for 
p ~ /i CT prohibits any further expansion of the scattering amplitude beyond 
(114) if the theory is to be applied for momenta < p < m^. Our phase 
shifts are given in Tab. 2. where for the reader's convenience we also included 
results of two modern potential models as given in [32,33]. The latter potentials 
reproduce the empirical phase shift rather well, where the Bonn potential 
shows a somewhat better x 2 as compared to the Nijm93 potential. We conclude 
that our chiral description of the ^o-phase shift does compete in quality with 
modern nucleon-nucleon potentials. 

It is instructive to discuss the singularity structure of our scattering amplitude 
(114). In Fig. 3 we show the scattering amplitude f^So^Q) evaluated below 
threshold with p = i q+e. The pseudo-bound state pole at q ~ —8 MeV and the 
anomalous thresholds at q = ± m 7r /2 are clearly seen. For unphysical negative 
q <C —8 MeV (2nd Rieman sheet) the amplitude is rather smooth whereas for 
positive q (1st Rieman sheet) the amplitude shows a rather complex structure 
anticipated in our discussion of section 2.3. The figure shows a well defined 
peak structure in the imaginary part of the amplitude and a corresponding 
zero of the real part at p ~ i 225 MeV, which was roughly predicted by the 
effective range theory (see (30)). From an analyticity point of view it is now 
obvious that a successful effective field theory description of the x So-partial 
wave amplitude requires the incorporation of this singularity. This explains the 
failure of the CE scheme (see section 3.3) which includes the pseudo-bound 
state pole but misses to describe the complete singularity structure seen in 
Fig. 3. 

In Fig. 4 we show how our UCE-Q 1 fit to the 1 S'o phase shift extrapolates 
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XPWA93 


25.3 
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-23.33 
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-23.59 


r( 1 5 )[fm] 


2.54 
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0.451 


0.451 
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0.6964 


0.7490 


6, [MeV] 


-558.9 


-672.5 


-724.1 


MMeV] 


354.1 


353.8 


353.5 


z [MeV] 


8.04 


-45.84 


-37.04 



Table 1 

Set of parameters obtained by fitting the 1 S'o(np) phase shift of PWA93. 
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E Uh [MeV] 


PWA93 


no pions 


UCE-Q 


UCE-Q 1 


Bonn 


Nijm93 


1 


62.07± 0.03 


62.06 


62.06 


62.06 


62.12 


62.07 


5 


63.63± 0.08 


63.86 


63.79 


63.73 


63.67 


63.62 


10 


59.96± 0.11 


60.03 


60.04 


60.02 


59.91 


59.94 


25 


50.90± 0.19 


50.46 


50.63 


50.75 


50.82 


50.86 


50 


40.54± 0.28 


39.83 


40.02 


40.20 


40.44 


40.38 


100 


26.78± 0.38 


26.43 


26.45 


26.53 


26.27 


26.17 


150 


16.94± 0.41 


17.08 


16.98 


16.96 


16.17 


15.71 


200 


8.94± 0.39 


9.36 


9.24 


9.15 


8.14 


7.07 


250 


1.96± 0.37 


2.35 


2.27 


2.19 


1.40 


-0.45 


300 


-4.46± 0.43 


-4.45 


-4.42 


-4.44 


-4.45 


-7.18 


350 


-10.59± 0.62 


-11.30 


-11.13 


-11.02 


-9.65 


-13.32 



Table 2 

We confront our result (112) and (119) with the PWA93 solution for the 1 So(np) 
partial wave. The x 2 of (112) with 12.9 and of (119) with 5.7 is significantly improved 
as compared to the x 2 of 25.3 which follows if pion effects are switched off (i.e. 
qa = 0,771,,- = 0). The set of parameters is given in Tab. 1. The x 2 are 37.0 for 
Nijm93 and 12.2 for the Bonn potential. Note that here we omitted the first point 
at E\ab = 1 MeV since its inclusion would somewhat artificially worsen the x 2 ■ 

to higher energies. Here we confront our result with Arndt's SP99 solution 
[33] . We observe that in fact our UCE scheme appears better compatible with 
the PWA93 solution of the Nijmegen group [32]. Note that the zeroth order 
(dashed line) and first order (solid line) results give very similar predictions 
for the phase shift for a given set of parameters confirming good convergence 
properties of our UCE scheme. The main effect of the correction terms sit 
at low momenta p < 200 MeV. Most dramatic is the modification at low 
momenta if the pion effects are switched off with g^ — as demonstrated with 
the dotted line in Fig. 4. Note that of course a set of renormalized parameters 
leads to a better description of the phase shift even with qa = as is shown in 
Tab. 2. Our full result (solid line) starts to deviate significantly from the SP99 
solution for momenta p > 450 somewhat above the pion production threshold 
at p ~ 370 MeV. 

We conclude our work with a short explorative discussion of the next terms 
in our expansion. At order Q 2 one needs to include the irreducible two-pion 
exchange diagrams in the interaction kernel K n studied in [8,9,34]: 



K. 
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Fig. 3. The 1 5o-partial wave amplitude /[i5 ](«9) = \fsT[i So ](iq)/(87r) of (114) 
evaluated with the UCE-Q 1 parameters as given in Tab. 1. 



(122) 



The scattering amplitude T n> s is to be evaluated up to order Q 2 applying the 
L-counting rule. This amounts at most a thrice iterated kernel K n ^s- Moreover 
1/m correction terms need to be included. Also the interaction K a s receives 
a contribution 





K a ,s = g*,s(s) J*C 



+ 



X'X'X-X 



(123) 



from a new chiral vertex 



£ = I g W A 4 (iV 75 Ctt 2 N l ) (n* C" 1 r 2 (fin) n) + h.c. 



(124) 
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Fig. 4. The 1 S'o-phase shift evaluated with the UCE parameters at order Q 1 as given 
in Tab. 1. The solid line is the full result (UCE-Q 1 ) obtained from the fit to the 
PWA93 solution. The dashed line is the zero order result of the UCE scheme and 
the dotted line follows from the UCE scheme with gA = (no pions). The dots 
represent the SP99 analysis of Arndt [33]. 

part of the chiral Lagrangian. The pion production parameter g^> is uncon- 
strained by chiral symmetry and may be estimated by meson-exchange phe- 
nomenology. We point out that the interaction K a s receives only 2-nucleon 
irreducible contributions which are also irreducible with respect to pion lines. 
This is a crucial element of our UCE scheme which only then leads to a renor- 
malized inverse coupling function g] i 1 (p,m 7T ) free of any pionic cut structure. 
Moreover the two-nucleon propagator G in (5,31) is corrected by the leading 
nucleon self energy loop diagrams: 



G- 



+ 



"W 



^*£W 



(125) 



The important new element of the order Q 2 calculation are the pion production 
cuts of loop functions. They give rise to a complex phase shift for p > fi n with 
the pion production scale ji 2 , = mm n + m 2 /4:. We estimate the result of the 
UCE-Q 2 calculation by adding all pion production contributions of order Q 2 to 
our UCE-Q 1 result (114). At the given order there are three different sources 
for pion production. First there is the pion production cut of the two-loop 
bubble diagram J^' R (p, z ) (see (68)) with 
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J S(P> z o) 



9 2 a\ 2 /Qj\ 



ml 



' (2tt)< 



'A + A? 



(p 2 



7T— production 



/4) 2 



In 



p l 



K-P 



In 



zo 



K 



+ iirQ(p-ii 1 c] 



(p 2 - A ) 
p 4 



2\2 



+ o (g 5 ) 



(126) 



We wish to emphasize that the chiral expansion of the pion production terms 
is not performed applying the heavy baryon technique since it would destroy 
the proper pion production threshold. Rather we expand the spectral density 
in powers of (p 2 — fi 2 ) n identifying fi w ~ Q as a further characteristic small 
scale. This procedure keeps the correct pion production threshold. Here we 
encounter a clear advantage of the relativistic chiral Lagrangian over its heavy 
baryon representation. 

We turn to the second pion production term. It results from the nucleon self 
energy correction of the one-loop bubble 



rW 



JkkiP, z o) 



2 \2 



9 a* 





(127) 

■n— prod. 



with its leading contribution proportional to the expression found for the first 
term (126) with Jsr(Pi z o) = ^J^r{Pi z o)- We point out that the dominant 
term at order Q 2 arises from the pion production vertex (124) which induces 
the loop function: 



J vl(P> z o) 



A 



-5 



9a 9 



« 



>c<*>e< 



+ 



X3<*>3< 



1 7r— production 



2 rn 4 n P 2 / V-/4) 2 b I /4 



(27r) 4 m\ v p 4 

(P 2 ~ A? 



Vl~P 2 



+ 



(*g + AY 
z Ip 2 



In 



/4 



+ Z7T 0(_p — y^Tr) 



Jf 



O Q 



A + z o 



128) 
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Technically the three pion production loop functions are incorporated easily 
by the formal replacement in (114): 

9r 1 ~> 9r 1 ~ 9 a A 2 4!fl(P. z o) - 9 2 a A 2 J^r(p, zo) 

-9R9Ag M X 5 4%(p,z ). (129) 

We point out that all loop functions J^ R (p, z ), J R (p, z ) ~ Q 4 and Jy R (p, z ) > 
Q 6 contribute to the renormalized inverse coupling g"^ 1 at chiral order Q A : we 
discover here the presence of the enhancement factor g R x ~ Q~ 2 which implies 
that at leading chiral order pion production close to threshold is driven by the 
pion production vertex (124): 



cw /^ ^ 9A(*rnj (p 2 - up (p 2 - pi) 2 , 

g\m\ (p 2 - p 2 a ) 2 (p 2 - p 2 w ) 2 ( x 

+ 2 WU? <W^ ^ p -^ (130) 



with: 



X A 9 M 9r 1 = ^-^- 2 +--- ■ (131) 

2 vr p 2 - pi 

In (130) we dropped some overall normalization induced by pionic vertex 
corrections for simplicity. Our result clearly demonstrates a subtle interplay 
of the pion production scale p n and \x a . Close to the chiral limit with p n <C \x a 
the second term in (130) is dominant for the threshold region with p ~ p n . 
For a realistic pion mass with p^ ~ p ff we find the first term in (130) most 
important since it is suppressed by only one power in p 2 — p 2 as compared to 
the two powers of the second term. 

In Fig. 5 we confront (130) with Arndt's SP99 analysis [33]. The inelasticity 
parameter p with 5s 5 = — In cos p as well as the real part of the phase shift 
are described well with ^ 2) = 0.751 fm, ^ 0) = 0.222 fm, £ CT = -783.5 MeV, 
p a = 350.4 MeV, z ~ -35.27 MeV, and ( w = 0.301 fm. Note that the set 
of parameters is rather close to the values given in Tab. 1 indicating once 
again the stability of our expansion. The pion production parameter (^ turns 
out natural since it implies g^ ~ —0.02. Fig. 5 shows also results of the 
perturbative expression (130) with ( n = 0.301 fm (dashed line) and ( n = f m 
(dashed-dotted line). This demonstrates that pion production is driven by the 
vertex (124) and furthermore that the production turns non-perturbative for 
momenta p > 400 MeV. 

The pion production cross section was studied close to threshold by a series 
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Fig. 5. The left figure shows the inelasticity parameter p with ^s 5 = —In cos p and 
the right figure the real part of the 1 5'o-phase shift. Our partial UCE-Q 2 calculation 
as described in the text (129) is given by the solid line. The dashed line is the result 
of the Bonn potential and the dots represent the SP99 analysis [33]. The dotted and 
dash-dotted lines follow from the perturbative expression (130) with ( w = 0.301 fm 
and C?r = fm respectively. 

of recent works [39-44] in the context of chiral perturbation theory. Since the 
total pion production cross section 



TV 



®tt— production 



2j9 2 



J2 (2J + l)(l-exp(-4 9<J Jl7 ) 

u 



(132) 



receives contributions also from higher partial waves even close to threshold 
[45], the comparison is not directly possible at this stage. Nevertheless it is 
amusing to observe that our production mechanism contradicts conventional 
chiral attempts to describe pion production insofar that in our scheme the 
production vertex (124) contributes at leading rather than at subleading order 
[39-44]. We therewith confirm phenomenological studies which demonstrate 
the importance of the a- meson exchange [46,47]. For a review of one-pion 
production in nucleon-nucleon collisions see [48]. 

We find our result for the x Sq partial wave encouraging and therefore the 
unitary chiral expansion scheme should be applied also to the spin triplet s- 
wave and higher partial waves. In particular we expect that the zero of the 3 S'i 
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and 3 Po-phase at £i a b — 300 MeV and £i a b ~ 200 MeV respectively is driven 
by a mechanism analogous to the mechanism presented for the 1 Sq phase. 
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A Partial wave representation of invariant amplitudes 



The on-shell nucleon-nucleon scattering amplitude is given in terms of five 
isospin zero amplitudes a\ (t,u) and five isospin one amplitudes a\ (u,t) as 
follows: 



T c (u,t)=V X>i 0) (M)7i + Pi J2af-\u,t)Ti 



i=0 



i=0 



(A.1) 



with the Dirac tensors 



Ti = (7 5 C- 1 xC7 5 ) , 
T 2 = ( 7a C- 1 xC 7 a ) , 
T 3 = (C- 1 xC) , 

T 4 = ( la C- 1 xC 7 j-| (a^C- 1 x Co**) , 
T 5 = (75 C- 1 x C75) + (757 Q C- 1 x C 75 7 a ) 
and the isospin projection operators: 



(A.2) 



V = \ (r 2 x r 2 ) , 
Vi = \ (t¥ 2 x T 2 f) . 

The tensor product 'x' in (A.2) is defined via the replacement rule 



(A.3) 



r c- 1 x cr - (u(pi)f c- 1 ^)) (V ( P2 ) cr u ( Pl )) (a.4) 

with the incoming momenta pi j2 and the outgoing momenta p[ 2 . The invari- 
ant amplitudes a) (u, t) can be uniquely reconstructed from the nucleon phase 
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shifts [35]. It is convenient first to express the amplitudes in terms of auxil- 
iary amplitudes j\ which have a more transparent expansion in partial wave 
amplitudes: 



m 2 m 2 m 2 u — t m u — t 

a l — 2 7l J3 — 74 7= —7 J 5 , 

S S S U + t J-S U + t 



m 

a 2 — 7= J5 , 

m 2 2 u — t m u — t / 2 \ 

«3 = 2^— h + m - — —-a / 4 + -7=7 -2 ( s + 4m j / 5 , 

m 2 „ m,\/~s 



a 4 — — Ja — - h , 

u + t u + t 



2 

m 



For the reader's convenience we also provide the inverse relations 



s u — t u + t 

/ 1 — 7, 2 fll ~~ O 2" ° 4 _ O 2" ° 5 ' 

2 m 2 2m^ 2m^ 



■U — t M + t M — £ 

72 = 2 — — a 2 + — — a 3 - — — 

u + t I m z 2 m A 

73 — 5~«5 , 



m 2 
. s u + t 

J 4 — — o «2 5~ 04 , 

m 2 m 2 

h = -^a 2 . (A.6) 

m 



The amplitudes fi are now easily expressed in terms of coupled 7n , / 2 2 , /12 



and uncoupled /j , / partial wave amplitudes 



/ 1 = E(2j + l)/o°' ) 4 ) W, 
i=o 

oo 

i=o 



+ (« + «)g w+ l)iff» + ^?) 



4 } (g) 

2^ 



/ 5 = -(* + «) £ (2j + 1) A2 ir^r ( A - 7 ) 



where d^, m (6) = (jm'\e~ l6Jy \jm). The sum over the total angular momentum 
j, in (A. 7) is constrained by proper selection rules. In (A. 7) we restricted the 
cm. scattering angle < 6 < ir for technical convenience. 

The coupled amplitudes /jy , f^i and fi 2 are expressed in terms of the ampli- 
tudes fi=j+i, fi=j-i with well defined angular momentum /,and the transition 
amplitude f± as follows: 



,o-)_ i+j , , j , , yi(i + i) , 



/22 - 2j + 1 ; ^' +1 + 2 j + 1 h=3 ~ l 2j + 1 J± ' 

A? = -^±S /,_, +1 + V^±^ /,_,_, + -J- /± . (A.8) 

j 12 2j + l J 2j + l J 2j + l 



The partial wave amplitudes are given in terms of the phase shifts. For the 
uncoupled amplitudes /q and /}•" one defines 



m z V M + tt v y 



with the partial wave phase shift 8j. The coupled partial wave amplitudes 



4vT / ~S 1 / / s } 2i5? 

I cos I 2 r ' ■ 

m J v u + t i 



^-S^tM^K^-i 



/f = ^/— sin ( 2 ef) e *(f=i-i+«f=i + i) (A.10) 

J m 2 \ju + t v J / v ; 

are expressed in terms of the phase shift $L- ±1 and the mixing angle e^ accord- 
ing to Stapp's convention [49]. Note that the spin singlet partial wave phase 
shift S^Sq) sits in fo. 
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B Evaluation of loops 



In this appendix we derive analytic results for the one, two and three-loop func- 
tion needed for the nucleon scattering amplitude at leading chiral orders. Since 
we apply a relativistic version of the chiral Lagrangian the non-relativistic ex- 
pansion in 1/m is performed at the level of the loop integrals. Technically 
this can be achieved by expanding the spectral density in powers of 1/m and 
then recovering the real part via a dispersion relation. Equivalently one may 
perform the energy integrations in the loop function first and then expand the 
integrand before working out the remaining three momentum integrals. When 
dealing with fermions a further complication arises due to the spin structure 
of fermions. However, any diagram involving complicated Dirac structure can 
be expressed in terms of scalar master loop integrals which naturally arise 
in scalar quantum field theory. This reduction can be achieved purely alge- 
braically relying only on covariance and the shift invariance of the integration 
measure dH (as implied for example by dimensional regularization). We there- 
fore exemplify our techniques first at hand of these master loop functions. 
Consider the scalar functions: 



hAs) = -i j -— S N (P - l) S N (K + I) , 



(27T) 

dH 

(27) 
hAs, t) = -i I -±L S v (l + \ Q) S v (l - \ Q) S N (P - I) S N (K + Z)(B.l) 



h,N(s) = i J -^- S n (l + l Q) S N (P - I) S N (K + J) , 



(27 

with Si(p) = l/(p 2 — mf + ie). For technical convenience we introduce the four 
vectors P, K and Q as follows: 



pi = K + \Q , p ' 1 = K -\Q 

P2 = P-\Q, p' 2 = P+\Q. 



(B.2) 



The master functions are easily calculated by means of dispersion techniques: 



4(m 2 +A 2 ) 4(m 2 +A 2 ) 

T , \ f ds' p 2 , N (s') f ds' P3,n(s') 



US' — S J 7T S' — S 

4 m 2 4 m 2 

4(m 2 +A 2 ) 

IaMs , t)= [ Hei^A (B .3) 



71 S ' — S 

4 m 2 
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with (P + K) 2 = s. The finite cutoff A < m is kept at this intermediate 
stage in order to mathematically justify the 1/m expansion and provide an 
instructive guide how to define a suitable regularization scheme. The cutoff A 
will ultimately be removed and may be replaced by an appropriate scale arising 
from dimensional regularization. We derive the s-channel spectral densities: 



m z 



1 /l 
1 1 



P2 '^ = 8^V4-V 



P3,n{S) 



16vrV4 



m^ — u 



t yJ—U — t 



In 



ml 



u — t 



mi 



P4,N(s,t) 



4tt v/=i y/Am 2 -u-t yjb(s,t) 
with u + t — 4 m 2 — s and 



arctanh 



t(u + 1) 



(B.4) 



b(s,t)=4:ml-4:ml(u + t)+t(u + t) . 



(B.5) 



In the algebraic reduction of the pion nucleon box diagram we will encounter 
further master loop functions: 



h,ir{t)=i 

hA t ) =i 



dH 
dH 



S N (l+p)S 7C (l), 



(2vr) 4 

dH 



S*{l + \Q)S«{l-\Q), 

-S n {l + \Q)S^l-\Q)S N (l + K) 



(B.6) 



We point out that the master loops of (B.6) are not independent of the master 
loops introduced in (B.l). They are interrelated by exact algebraic identities 
required for the cancellation of kinematical singularities: 



h{m 2 ) = I 2 , N (4m 2 )-mlh, N (4m 2 ) , 
hAt) = hA^m 2 -t)- (ml - \t) h, N (4m 2 -t,t), 

2 =/ 3 ,7v(4m 2 -t) 
\Um\- t) 7 4) jv(4 m 2 - t, t) + 2 (2 m\ - t) I' 3N (A m 2 - t) 



2m 2 n ~t) §- s h, N (s,t) 



f (ml h, N (4m 2 -t)- I 2 , N ^ m -t)+h\ 



m 



Am 2 -t 



ml- It) hAt) + hAt) ~ h(m 2 )) . 



(B.7) 
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The identities (B.7) can be verified for example in dimensional regularization 8 . 
Any consistent regularization procedure must respect (B.7) order by order in 
the 1/m-expansion. Once we specify our regularization for /^(s) in terms of 
the cutoff A the regularized loops hip 2 ), h^if) an d h^{t) must be compatible 
with (B.7). 

One may conjecture that the triangle loops l2,n{t) and l3,n(t) could be evalu- 
ated by means of a dispersion integral: 



A[ml+K 2 ) 4( m 2+A 2 ) 

with the t-channel spectral densities: 



e(t-4ml) [i ^ 

P2,w(*) = — 




t ' 



p^(t) = — , V wJ arctan y/Am*-t V - (B 9) 

8tt ^/t(4m 2 -t) V t-2m2 ^ 

for t > 4mJ. The pion bubble h^if) is rendered finite by one subtraction. 
Performing the dispersion integral for l2,n{t) in (B.8) one finds: 



hA^ m l) = 3~T arc si n h 

87H \m 




Um 2 —t / £ 

/ 2 ..(f) = W0)- — 1 \/— =— arctan ^-j— - ] j (B.10) 

where we applied the large A-cutoff limit in the subtracted loop. We point out, 
however, that the cutoff A in l2,n{t) and l3,n(t) must not be identified with A. 
The 1/m expansion of the loop integral l3,n(t) is subtle. The regularization of 
h,n(t) by means of the cutoff A <C m is incompatible with the 1/m expansion. 
This follows from the fact that the 1/m expansion of the t-channel spectral 
density leads, at subleading orders, to singularities in the dispersion integral 
at t' = 4mJ whereas h, n (t) given via (B.7) leads to well defined moments in 
the 1/m expansion. 



8 More conveniently the identities (B.7) can be derived algebraically if one requires 
the reduced loop function bis(s,t) of the pion-nucleon box diagram not to be sin- 
gular at u = and u + t = 0. See appendix B.6. 
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At the one-loop level we complete the collection of master loops with the 
pionic and nucleonic tadpole diagrams I\^ and I\n'- 



4(m 2 +A 2 ) 



4 m 




(2tt) 4 P -ml + ie 



A J ml + A 2 - ml arcsinh ( ) ) (B.ll) 



. 7r 2 \ " " \m 



where consistency of the cutoff regularization relates I\^ with ^s I 2 ,n{s) and 
7i ]7r with $5 l2,-n{t). It is crucial to observe that the proper A dependence of the 
pionic tadpole ii ;7r follows from a further exact algebraic identity 



h,* - h,N = 6m 2 ml [ml — I^ N (4m 2 ) - — J 2 ,Af(4m 2 ) 

\ as as J 

+ (2m 2 ml + ml) I 3 , N (4m 2 ) + (m 2 - m 2 n ) / 2 , ff (4m 2 ) (B.12) 

which may be verified in dimensional regularization. Again we insist that 
(B.12) must be respected by the 1/m-expansion. In particular (B.12) implies 
a non-trivial matching of the s-channel cutoff parameters A and the t-channel 
cutoff A. 

We point out that all t-channel loops Is y7T (t), l2,n(t) and 7i ;7r can in fact be 
defined via (B.7) and (B.12) in terms of the s-channel loops I 4jN (s,t), ^^(s) 
and I2,n(s)- It is, however, convenient to simplify our regularization scheme 
and derive more explicit expressions for the loop functions in order to facilitate 
the renormalization procedure. First we assume that the diagrams are already 
partly 1/m renormalized via the construction (46). This implies that the factor 
l/v^s 7 in the spectral functions p^n{s') can be pulled out of the dispersion 
integral (B.3). This procedure moves a large class of power divergencies of the 
form (A/m) n into the interaction kernel where they are absorbed into the local 
interaction terms. Next we systematically drop nucleonic tadpole diagrams 
since they also involve multiple powers of the cutoff A™ with n > 3. Such 
terms are ill-defined and are ultimately canceled from contributions induced by 
the properly constructed cutoff dependence of the path-integral measure [27]. 
Finally we apply the large cutoff limit for all finite parts. Our regularization 
scheme leads now to at most linear divergent one-loop functions. 



63 



In more formal terms our regularization scheme is defined by writing the one- 
loop functions I 2 ,n(s), h,N(s) and I 4jN (s,t) as follows: 



I 2N (s) = - -= l-A + ip) +A/ 2 ,jv(s) 



I 3tN (s) = —^ - hi (l - 2 i-?-) + AI 3 , N (s) , 
16 7r y/s p V m w J 

Ian(s, t) = ■= = i arctanh 




1 ( y/-tb(8,t) y/=£=t (, 

-\ — arctan | - — - — 1 — 

Ami rriir \ 



yj-tb(s,t) y/-u-t ( 
Am\ m n \ 

AI^ N (s,t). (B.13) 



1 

+ - arctan I - — - — h 

2 \ Ami m 




where 



4(m 2 +A 2 ) 

AT /x f ds"SI 2 , N (s' 

M 2 ,n{s) = - / 

J IT S - 



- Vss' 

4 m 2 



4(m 2 +A 2 ) 

AI 3tN (s) = - 



ds' $s I 3>N (s') 

7T s + y r ss' 

4ra 2 
4(m 2 +A 2 ) 

Ah N (s,t) = - f d A^M0. (B .14) 

4 m 2 

Note that for example the separation of the term A/ 2 ,at(s) from I2,n(s) in 
(B.13) is induced uniquely by our construction (45,46). The divergent structure 
A/ 2 ,at(s) can systematically be absorbed into local counter terms in accordance 
with (46). This follows from the representation in (B.14) which shows that 
AI 2 ,n{s) can be expanded in a Taylor series around s = Am 2 . In our simplified 
regularization scheme such terms are from now on simply dropped. The loop 
functions I 3 ,n(s) and I ijN (s,t) require special attention since AI 3jN (s) leads 
to chiral log terms and AI^^(s,t) depends on t in a non-polynomial fashion. 
It is instructive to consider first AI 3jN (s) and AI ijN (s,t) as defined via (B.14) 
in terms of the finite cutoff A and inspect the point s = Am 2 , t — 0: 



* T / 9n a l i /2A\ ir m w ,^/A 2> 

AI 3:N 4 m 2 ) = - 3 1 - In — - - -1 + O _ 
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A/4,7v(4m 2 ,0) 



A 



64 7T 2 m 3 ml 



1 + 



A 2 



TO 



(B.15) 



where A < to. Note that A/3 i Ar(4m 2 ) must not be considered as a counter 
term since it is non-analytic in the squared pion mass. Also the term A/m 2 
cannot be absorbed into local interaction terms of the chiral Lagrangian. The 
result (B.15) signals an infrared singularity induced by the 1/m expansion 
and asks for a more careful regularization procedure. In fact the solution to 
this puzzle follows simply if one evaluates A/4 ; jv(4to 2 , 0) for an infinite cutoff 
A > to: 



A/p(4m 2 



1 



lQir 2 m 2 



In 



/;/ 



m„ 



+ 71' 



m 



m„ 



-1 



4 m 2 
1 + ln 



— 1 arctan 



4 m 2 

\ m l 



- 1 



/;/ 



A/ 4 ,iv(4m 2 ,0) 



16 7r 2 m 2 \ m„ 

-1 fir 2m 

16 n 2 mm 3 I 2 

- 1 U 

32 it 2 m 2 m 2 \ 8 to 6 m' 2 



8 m 12 m? 



m': 



m° 



arctan 



4 m 2 



mt 



Am 2 

\ m l 



7r m-n 1 m 2 „„ ( mS\\ .„ 



This is a surprising result since it demonstrates that the chiral log term ln(m^) 
in (B.16) results in part from modes probing large nucleon virtuality with 
A > to. 

In order to perform the 1/m expansion of AI^n(s, t) and derive explicit results 
it is convenient to perform the angle average in (8) after the s'-integration in 
(B.14). Now with A = oo in (B.14) we obtain: 



A/ 4 , W (4m 2 ,{) 



1/2 

-1 f dx 



■ 7T^ 



7T 



arctan 



y-1 +a(x,t)J 



, , K y2y/a(x,t) 

a{x,t) = ^j (4+^-(4x 2 -l) j . 
m z n \ m^ v 'J 



yj—l + a(x,t) 



(B.17) 



The 1/m expansion of AI^^(s, t) follows upon expanding the integrand of 
(B.17) in powers of 1/a: 
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/arctanh 

L ' — ^ — tyyt ^ - ■ J l l • ■ ■ ■ j- ir 



iir* m^ 



arctann v / I ^^ ^ x g? f fe.Y' 

V^4m2-tv^t 8m \.m\-t) W7 



where £ < 4mJ. Note that the 1/m expansion of AI^N(s,t) stops converging 
at £ = 4mJ. 

We turn to the master loop function J 3j7r (£). It is instructive to provide the 
loop including next-to-leading-order contributions 



T .. 1 / t \ 1 /V-M m *- 2mi-t 

h*(t) = 1 + n -^= arctan -— £ 

^ w 16 7rm V 8m 2 / x /=t \2mJ 1287rm 3 4m 2 -£ 



1 2m 2 -£ / =T 

H — - . arctanh 



16 7T 2 m 2 / t(4 m 2 _ t ) \Aml-t 



+ A/ 3 , J v(4m 2 ) + o(-l) (B.19) 

which follow with (B.7,B.13,B.18). The result (B.19) is in agreement with the 
leading expression derived in the heavy-baryon framework [36,37] and with 
(B.6,B.13). The subleading term of order 1/m 2 agrees with a corresponding 
expression given in [38]. Our result (B.19) confirms that the triangle loop 
h,n(t) and the scalar pion nucleon box function I^^(s,t) show a singular 
branch point at t = 4 m 2 if expanded in 1/m. A systematic resummation 
leading to a smooth behavior at t = 4mJ is proposed in [38]. 

The function AI^^s, i) can alternatively be deduced from (B.7). Imposing 
the functional form for I2,n(s), h,N( s ) and I^ N (s,t) as given in (B.13) with 
Ai2,jv(s) = we derive: 



^MH^MiW^w^a). (B , 0) 

2 m 2 (4 m 2 , — t) \m 4 / 

Our results (B.20) and (B.16) determine the divergent subtraction h^ify in 
terms of i2,7v(4m 2 ). This is a surprising result since it shows that the 1/m 
expansion converts the log-divergence in I2,n(s) and l2,n{t) into a power 
divergence. In fact with (B.12) and (B.7) we find that also the pionic tadpole 
diagram 7i ;7r and hijn 2 ) must not be dropped and assigned a linear divergence: 



h,« = (m 2 -\ml) J 2 ,^(4m 2 ) + 4m 2 m 2 A/ 3 ^(4m 2 ) + O (^) , 

/ 2 ,.(0) = i/ 2 ^(4m 2 ) -2m 2 A/ 3 ^(4m 2 ) - -L - _^!L_ + o ( ^) , 
2 8 it 2 64:7im \m 2 I 
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I 2 (m 2 ) = I 2 , N (4m 2 )-mlAI 3 , N (4m 2 ) - -p- . (B.21) 

We emphasize that our result (B.21) follows from the convenient prescription 
to systematically drop the divergent terms i^jv and A/ 2 ,jv(s)- This defines our 
simplified regularization procedure. The regularization of the remaining loop 
functions is then implied by the algebraic constraints (B.7,B.12). The latter 
follow for example in dimensional regularization and must be satisfied by any 
consistent regularization scheme. To this extend it is not surprising to find 
that for example the pionic tadpole in (B.21) acquires an explicit dependence 
on the nucleon mass. The only divergent structure left at the one-loop level is 
/2,7V (4 m 2 ). It may be written in terms of the cutoff parameter A (see (B.13)) 
or equivalently in terms of a scale resulting from dimensional regularization. 
Consistent power counting rules arise if we count /2,7v(4m 2 ) ~ A/m ~ Q. We 
conclude that the scalar nucleon bubble is in fact the only master loop where 
we have to strictly insist on A < m since for A ^> m one finds I 2) n ~ In A ~ Q°. 
All finite loop functions h^it), h,N(s) and I^N(s,t) must be evaluated with 
A ^> m in order to generate the proper chiral log content and avoid infrared 
singularities (see also [50]). 



B. 1 Static pion limit 



Here we would like to confront our relativistic scheme with the potential ap- 
proach. Inherent of the potential approach is the use of the static pion ex- 
change. The pion propagator is expanded 



where it is commonly argued that terms involving the energy transfer Z 2 , are 
suppressed by 1/m. Mathematically it is not quite obvious that such an ex- 
pansion is justified and in accordance with covariance since the expansion is 
highly divergent due to multiple powers of the energy transfer / . We evaluate 
the master pion nucleon box function I^^(u, t) in the limit of static pions 



j-(static) / -> -,/\ _ 1 f d I 1 1 1 

A ' N [P,P ~ 4^ J (2tt) 3 (f+ py + m 2 ( f + p y + m 2 p - p 2 - % e 

where we chose the center of mass frame with p 2 = p' 2 and p ■ f? = p 2 cos 8. 
Note that we dropped terms suppressed by 1/m 3 . We compare the static pion 
result (B.23) with the proper non-relativistic expanded master loop function 
of (B.13) and find agreement. 
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B.2 One-loop bubble 
Start with the one-loop bubble function 

dH 



J(s) = -itr J—-C l5 S F (P + l) l5 C~ 1 S t F (K-l) , 



A 2 



1 f dp 2 pj(p) mp I p 2 

IT J p z — p z — %e Z7T V w 



^.A) = r/S^T,P,M = ^Mi + S (B-23) 







with the nucleon propagator Sf(p) = l/(jf> — m + ie). The loop function may 
be regularized by a cutoff A introduced in the s-channel dispersion relation 
with s = (P + K) 2 = 4 (m 2 + p 2 ). Equivalently the loop may be expressed in 
terms of the master functions: 



J(s)=2s/ 2i7V ( S )+4/ 1)7V . (B.24) 

The expressions (B.23) and (B.24) agree if the cutoff regularized tadpole is 
given by (B.ll). Our simplified regularization procedure leads to a linear di- 
vergent one-loop function J(s, A): 



ml P 2 /"2 

— Jl + ^2 - 

2ir V m A \7T 



J(s,A) = -\ll + ^[-A + ip) . (B.25) 



At given order in the \jm expansion the renormalized loop function is intro- 
duced by the minimal number of subtractions at p = — i z. At leading chiral 
orders this leads to: 



Jr(p,z) = —- (-z + ip 

Ztt v 

Note that the leading term is renormalized via Jn(—iz,z) = whereas the 
subleading term requires a further renormalization condition. 



B.3 One-loop vertex 

The one-loop vertex function is introduced with 
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i r dH 



■S K (l + ±Q)S F (K-l)f>y 5 U+±<8), 

u t (p 2 )V^(p 2 ,p 1 )u(p 1 ) = V n (s)u t (p 2 )C^T 2 T t u( Pl ) (B.27) 

where u(p\ )2 ) are the incoming nucleon spinors. The loop function is expressed 
in terms of the master functions 

VM = 4 m 2 (ml I 3:N (s) - I 2 , N (s)) + 7 1;7r - 2 I hN + 2 m 2 n I 2 (m 2 ) . (B.28) 

According to our renormalization procedure we first transform to the on-shell 
irreducible interaction kernel. This leads to the pion-exchange with pseudo- 
scalar coupling which is known to be on-shell equivalent to the pseudo-vector 
coupling used in (B.27). The vertex evaluated with pseudo-scalar coupling 
agrees with the pseudo-vector result (B.28) up to the divergent terms ii ;7r — 
2ii,7v + 2mll 2 (m 2 ). The latter terms are included in the local part of the 
on-shell irreducible kernel K = K n + K a according to (44): 



Ka = Us) + ^fi- fa,* " 2 h,N + 2 ml I 2 (m 2 ))) T 1 V 1 + ---. (B.29) 

with the bare coupling function g(s) of (14). Similarly the AI 3y N(s) and 
A-/2,iv(<s) terms in (B.28) may be absorbed into K = K n + K a according 
to (46): 



K«={g{s) + 4^ (h,,-2I 1 , N + 2mll 2 (m 2 : 



2fl 



,2 



+ ^^4m 2 (mlAI 3 , N (s)-AI 2 , N (s))\T 1 V 1 + --- . (B.30) 



Note that in our regularization prescription the nucleonic tadpole i^jv is 
dropped whereas the pionic tadpole (see (B.21)) contributes at chiral order 
Q. The renormalized vertex function K-,_r(s) follows upon a proper loop sub- 
traction according to (41,63) with 

VM = V K>R (s) - c(s) J(s) + Ji )W - 2I 1>N + 2mll 2 (m 2 ) , 

t/ f\ m 2 m n m n ( . 2p\ 1 m 2 

V k ,r(s) = - , 5 1 In 1 - i , c{s) = - . (B.31) 

8 7rv m +P V V m ?r/ 2 m A + p z 

Note that the loop function K-,.r( s ) is most easily derived by first evalu- 
ating its imaginary part. The real part is then constructed with ease since 



69 



\/m 2 + p' 2 K-,h(4?71 2 + 4p 2 ) is an analytic function in the upper complex half 
p-plane. 



B.Jf. Two-loop bubble 

The two-loop bubble integral is 



Us)=~tr J -^TjT^C-'S^P + l)V«\P + l,K -l)S F {K -I) 

(B.32) 

defined in terms of the full off-shell one-loop vertex function (B.27). Here 
we implicitly assume that the vertex (B.27) is evaluated with pseudo-scalar 
pion nucleon coupling according to (44). We evaluate J v {s) by constructing 
its imaginary part and then recovering its real part by means of dispersion 
relations. The loop is written as a sum of three terms 



Ms) = Us) - c{s) J 2 (s) + 4*>(s) (B.33) 

where J^(s) is the unique contribution due to the pion production cut. The 
square of the one-loop bubble J(s) in (B.33) takes care of the overlapping 
divergence structure of the two- loop function J^s). 

We first evaluate the loop-subtracted bubble function J n (s) by imposing a 
dispersion relation on s J n (s). Here we pull out the factor 1/s' in the dispersion 
relation in order to comply with our regularization scheme and (46). With c(s) 
as given in (B.31) the appropriate spectral density is readily derived 



%.Us) = 2pj(s)$lV n , R (s) (B.34) 

in terms of the renormalized vertex function K-^(s) of (B.31). The following 
representation can now be established 



J M - J.(0) = i ^ V<» = -2 ^ In (l - i m (B.35) 

n ' (An) 1 \ m w J 

where y/sV^^s) = mV^ #(p). Note that (B.35) shows that all \/m correction 
terms for J n (p) — .7^(0) cancel identically. The identity (B.35) follows from the 
Kramers-Kronig relations: 
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'MSoo-tfiotf """""^W) 



— oo 

7T 



7T p p — p — le 

Up) - Xr(o)) - <v$(o) , 



J -2 2 

J.(P) " J-(0) = / ^- ^ _ 2 T . K ^i(p) • (B.36) 

J 7r z p z p z — p z — ie 



This is an amazingly simple result. The two-loop bubble function is given in 
terms of the one-loop vertex function \Qr(p). Note that the loop function 
Jir{p) shows a log- divergence. In analogy to (B.23) the subtraction constant 
Jtt(O) may be reconstructed by a cutoff regularized dispersion relation: 



1 df jm p\ „.,$),., 



J.(0) = m»/^ ?+* RV«((«) 

J it z \ v ml 



7i^ \p m / 

The renormalized loop follows now with J ff) jj(p, z) = Jn(p)—Jn(—i z) according 
to (42). 

We turn to the evaluation of the pion production term J^(s). This contribu- 
tion is conceptually interesting since it involves a further small scale 



o 

& = mm n + -z- (B.38) 

characterizing the the pion production process at p > fi n . The spectral density 
is systematically expanded in powers of (p 2 — /i 2 )™ 



pg(p) = mi 6(p - /i.) ^ff + ^ (Q 5 ) (B.39) 

where we count p ~ fi n ~ Q. We point out that our expansion of the pion 
production spectral function is in full analogy to the non relativistic expansion 
of the two nucleon unitarity cut where one expands in powers of the relative 
momentum p also around the threshold value p = 0. The real part is derived 
from the dispersion relation: 



(2vr) 4 \ p 2 p 



2 



4'lM = ^ -^f + ^^m '" ;l 



^2-p 2 
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/4 ^ + A? 1n 

z 2 z* 



/4 



/4 + 2 2 



+ ^7re(p-/x 7r ) (p2 ^ 4 /i ' )2 ")+0(Q 5 ) (B.40) 

properly subtracted at p 2 = —z 2 . Note that the chiral order J^' R {p,z) ~ Q 4 
is deduced by applying simply Weinberg's counting rules. In other words due 
to the absence of an intermediate nucleon pair state there is no enhancement 
factor for the pion production spectral density. 



B.5 Pion production loops 

We consider the self energy correction to the one-loop bubble diagram 

Ms) = -itr J AL C-S S F (l + \w) E (l + \w) S F (/ + \w) 

l,CS F (-l + \w) 
-^ j -^ C ~ l ^S F (l + \w) lb CS F (-l + \w) 

E* (-/ + \w) S F (-/ + \w) (B.41) 



with w 2 = s and the self energy 



E(p) = E^-3^A^|— — l5 fS F {p + /) 75/^(0 +c.t. (B.42) 

receiving contributions from the pion nucleon sigma term E^jv — 45 MeV and 
suitable counter terms. We also consider the two- loop function, Jy(s), induced 
by the pion production vertex in (124): 



Ms) = -2 i tr J -^L c- 1 v(l + \w) S F (l + \w) 75 C S F (-1 + \w 
-2itr jALc-^Sp (l + \w) V (-1 + \w) CS F (-1 + \w) , 

V(p) = -tJ-^yJS F (p + l) l5 fS n (l) . (B.43) 

with the vertex function V(p) = 70^ (p)7o- Note that in (B.41, B.43) we 
suppress the isospin space but include isospin factors effectively. 
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The two- loop functions Jt,(s) and Jv{s) have two types of s-channel cuts. Cuts 
which do not separate the pion line simply give rise to a renormalization of 
the nucleon mass in the one-loop bubble J(s) in case of Js(s) and a renormal- 
ization of the bare coupling in case of Jy(s). They are effectively taken into 
account by using the physical nucleon mass in J(s) and a renormalized bare 
coupling. The more interesting cuts with the pion line separated reflect the 
pion production process. We derive the leading contribution: 



P { jkp) = 3 < 6(p - /i.) ( £L fK + O (Q l 



2(2tt) 3 p 4 

and conclude Jj?r(p, z) — 3 J„ R (p, z) + (Q 5 ) (see (B.40)). Again we confirm 
that Weinberg's counting rules can be applied for pion production cuts. 



B.6 Pion nucleon box 

Consider the pion nucleon box diagram 
alH 



B = ~ % J JW SAl ~* Q) Tl1 " ( /+ ^) S ^ P ~ ^ 5 ("/+ 50) 

®r il5 (/+§$) S F {K + l)T jlb (-l+\Q) S 7r (l + lQ)(BA5) 

where the tensor product '(g/ is defined via the replacement rule 

f <g> r -► (u( P ' 2 ) r u( P2 )) (tz(pi) r u( Pl )) (b.46) 

with the incoming momenta p\^, and the outgoing momenta p[ 2 of the nucelons. 
On-shell projection simplifies the diagram considerably. The isospin one com- 
ponent, B^\ follows with: 

B (i=D = x ^ (^ lv) + ^ ® X » + X»® lll )+X (B.47) 

where 

X^ = lQmH J ALl^S^l - \Q) S N (P - I) S N (K + I) S«(l + \Q) , 



(2k) 
dH 

(2k) 



X, = Am 2 ij ALl IJt (j + 2m)s^l-\Q)S N {K + l)S 1x {l + \Q) 
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X^Am 2 ij— yj^l- 2m)S„(l-\Q)S N {P -l)S v (l + \Q) , 

x = ~ l S 4^r ( 2 m - /) ® ( 2 m + /) s *( l ~ \ Q) s *( l + \ Q) ■ ( B - 48 ) 



It is useful to introduce auxiliary loop functions Ji(t) and 6j(s, t) which exhaust 



the structure of X^, X^, X v and X: 



X, u = b 1 (s,t) 9 ^ + b -^l(P + Ky(P + K 

X, = ~\ J?\t) 7, - (| 4*\t) + 2 Jj*>(t)) ^ , 

*„ = -| ^(«) 7, - (| ■/J'ty) + 2 JW(*)) ^ , 

^=5 jr o* ) (*) (7^7^-2^^) (1®1) . (B.49) 

Here we drop terms proportional to Q^ or Q u since they would not contribute 
to the on-shell scattering matrix. Note that transforming to the on-shell irre- 
ducible one-pion exchange with pseudo-scalar coupling 'moves' all Jj func- 
tions into the on-shell irreducible kernel. The invariant functions bi(s, t) follow 
easily after applying the substitution rule 



11 ^a ( l 2_(l^l_(±Wl_(±iK-P)r\ (B5Q) 

W,W u ( (l-Q) 2 (l-W) 2 (l-(K-P)f \ 
W 2 [ Q 2 W 2 (K-P) 2 J 

| (K-P),(K-P) u f (l-Q) 2 | (l-W) 2 , n (l-(K-P)) 2 ]2 \ 



(K-P) 2 \ Q 2 W 2 (K-P) 2 J 

in (B.49) with W = P + K which is justified by covariance. The reduced box 
functions, bi(u, t), can now be expressed in terms of the master loop functions 
systematically: 



MM) = 16 ro 4 | ym \ *' -ml + \\ I 4 , N (s,t) 

+ l 6m 4 ( l <^\\ J (s) + 16m 4<-li / (t) 

\ 2 u I u 
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(m 2 



t^ 



b 3 (s,t) = 16— \ml-{-2 ^ " " | I 4 , N (s,t) 
u \ u 



m 6 / 1 ml _ ml — \ 

'3,7V 



i6— [;-^2L T -2^^^ / S J S ) 



■u \ 2 -u + £ -u 

2- 



m v m 2 — 4 T „ s . m 



-32— ^2— 1/^(4) - 16-^— (l 2 , N (s) -I 2[ m 

u u u[u + t) v 

-16 12- I^Lzij 3 ( t )_i 6 _J^(/ 2jW (t) - J 2 (m 2 )) (B.51) 

7/. .S 4- 11. II. ft -+- 11, 1 \ / 



■U S + M " w(s + Ml 



with m = 4 m 2 — s — t. The box function 62(5, t) can in fact be written in terms 
of bi(s,t) and j\ (£) as follows: 



2 2 

6 2 (s,t) = -— 6i(5,t) + 2 — Jf } (t). (B.52) 



It remains the evaluation of the reducible loops Jj (t) . Applying similar meth- 
ods we derive: 



4"\t) = lh,,-(ml- t I )hAt), 

j«"( i )=4^('M_il_ m j + '| /3 , w 

s + u 4 



Am 2 
+ 



s + u 



(ml - |) (/ 2 ,.(t) - / 2 (m 2 )) + 2m 2 / 2 (m 2 ) , 



2 1 ^ 



s + -ul 4 S + M 

16 m 4 / 2 i\ 

16m 4 / m 2 -±£ 1 m 2 \ _ . 2 . 8 m 2 /r r \ 

+ ^_ 3 2 -- + — M J 2 (m 2 ) + — — (/ 1>w -7 1;JV ) , 

Ji 7r) (t) = -2m 2 / 2 , w (t) 

jW( t ) = 8 _?£_ ( m 2 _ 1 t ) 7 (t) + 8 ^L h (*) _ / 2 ( m 2)) (B.53) 

with s + w = 4 m 2 — t. Finally it is convenient to decompose the on-shell 
projected pion nucleon box into our Dirac tensors T i} with simple exchange 
properties 
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B(p' 2 ,p 2 ;p' 1 ,p 1 )=V £ BP(u,t)T i + P 1 ]T B[ 1] (u,t)T t . (B.54) 

i=l i=l 

The exchange box diagram is now deduced by properly exchanging u and t in 
the reduced functions Bi(u, t) (see (56)). Note that the isospin zero amplitudes 
are B\ (u,t) = %B\ (u,t). Here we provide the necessary on-shell tensor 
identities (see [24]): 



(fi- m )®( m -f>) = {±t-lu)T 3 + \(u + t)T A + \(u + t)T 5 , 
$ <g> m + m ®f = (m 2 - \(u - *)) 71 + (m 2 - \(u + tj) T 2 

- (| m 2 + \ [u - tj) T 3 - (i m 2 - \ (u + tj) T 4 

+ |m 2 T 5 , 
$ <g,f = (I w 2 - \ (u - tj) 71 + Q m 2 - \ (u + tj) T 2 
-(\m 2 -\u)Tz-\m 2 T A 

+ (|m 2 -i( M + t))T 5 , 
ig + m) <g> (m +f) = (2 m 2 - \ (u - tj) T 1 +(2m 2 -i( U + tj) T 2 

- (m 2 - \{u + tj) (T 3 + T 4 - T 5 ) . (B.55) 

We collect the leading terms of chiral order one: 



B[ 1 \u,t) = -2m 2 (ml - f) I^{t) - * fe^t) 



2/ ■Vv'V 4m 2 

+ 2m 2 / 2 (m 2 )-ii/ 1 ^ + --- , 
S?' («,*) = "2 ™ 2 (< - I) 4- (t) + 2 m 2 J 2 (m 2 ) - g J 1)W + • • • , 
B^\u,t) = -\b 1 {s,t) + t T ^h{s 1 t) + --- , 
Bj l \u,t) = {b 1 (s,t) + f£b 3 (s,t) + ■ ■ ■ , 
^ 1) ( M ,t) = i6 1 (s,t) + |±|63(s,t) + --- (B.56) 

where the dots represent further terms with chiral power ~ Q 2 not relevant 
in this work. 

It is instructive to consider also the pion-nucleon box diagram defined with 
a pseudo-scalar pion-nucleon coupling vertex (PS). In our calculation the 
pseudo-scalar box follows formally from the pseudo-vector box diagram de- 
fined in (B.45) by dropping the X^,X M and X terms in (B.47). We observe 
that the difference of pseudo-vector and pseudo-scalar box diagram carries 
leading chiral order Q°: 



76 



BP(u,t) - B[^ PS \u,t) = -Am 2 (ml - f ) / 3>w (t) - 2m 2 / v (i) 



S^V,*) - B£' PS) (u,t) = -Am 2 (ml - §) / 3>w (t) - 2m 2 / v (i) 



+ 4m 2 / 2 (m 2 )-i|/ 1 ^ + 0(g 2 ) . (B.57) 

This implies that the on-shell irreducible kernel K = X^f=i 9i( u i ^ ^1 °f (44) 
receives a non-local contribution. For example: 



gi(u,t) = ± g\ A 4 (B? (u, t) - B^ PS) («, t)) 

+ ^lA 4 (5i 1) (t, M )-i?i 1 ' P5) (t, M )). (B.58) 



The non-local terms of (B.58) are then included in K n where K = K n + 
K ff . The pionic tadpole in (B.57) is included in K ff and absorbed into the 
coupling function g(s). Since both K n and K n are treated in perturbation 
theory our procedure (44) defines a systematic procedure to absorb divergent 
terms, which are not associated with the s-channel unitartity cut, into the 
local bare coupling function K a . At sub leading order it leads to the inclusion 
of the pionic tadpole I\^ into g(s). 

We close our essay on the pion-nucleon box diagram with the relevant combi- 
nation for s-wave scattering in the spin singlet channel: 



rp I +\ (1)/ +\ U ~t (1)/ + \ U + t (1)/ + \ 

s s 



2g\ A 4 m\m 2 h, N (s,t) - Ag\ A 4 mlm 2 I 3 , N (s) 

1 

2 



+ 2 g\ A 4 m 2 7 2;JV (s) + m 2 O ( -^ ) + (u <-► t) (B.59) 



where we now dropped the pionic tadpole contribution. 

B.7 Two-loop vertex 

The two-loop vertex function is 

Gin /„ ,_ „ _ _ , _xt. 



v«(p 2 ,p 1 ) = ^/^ I (^(p + /)f 75 (/+i^)) t vW(p + /,^-/) 
-ss + \q)s f (k -i)T lb y+\(%) , 

u\p 2 ) VW(P2, Pl ) w(p!) = V^( 5 ) </(p 2 ) ((7 75 ^ r 2 r,) ufo) . (B.60) 
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We evaluate V n7T (s) by imposing a dispersion relation on the renormalized 
vertex sV W7T:R (s) with 

VUs) = V n7TtR (s) - c(s) J(s) V v>R (s) - c(s) (Us) + jW(s)) 

+ c 2 (s)J 2 (s) + VJ:] R (s). (B.61) 

Our renormalization via loop subtraction keeps track of the overlapping diver- 
gency structure of the vertex function V W7T (s). The unique contribution from 
the one-pion production cut is denoted by \Q*' R (s). With c(s) given in (B.31) 
the renormalized loop functions V W7r: ji(s) and V^ R (s) are finite. Note that 
the pion production contribution \Q„ R (s) is renormalized by the proper sub- 
traction of the pion production loop Jjf\s). From (B.40) we expect V^. R (s) 
to carry chiral order Q 4 . The spectral function of V nn , R (s) follows in terms 
of the renormalized one-loop vertex V^^s) and an effective pion-nucleon box 
B^Js): 



9V^r(s)= 3^(s) ) (UV^Ris)) + 3J( S ) WB 



dx 



/d x 
— h,N( s 't) 

-1 



+ O (Q 2 ) (B.62) 

t=-2p 2 (l-x) 



with s = 4 (m 2 +p 2 ). Note that the function B niT {s) coincides with the s-wave 
projected pion-nucleon box diagram R nn (p) (see (80)) as should be expected. 
The renormalized loop is decomposed into its chiral moments: 



W) = E *&(P) (B-63) 

n=2 

where V^ R (p) carries chiral power n. The real part of the two- loop function 

(2) 

V^ n R (p) can be easily derived without performing further integrations. Recall 
that V^n r{p) satisfies a dispersion relation in p: 



VS M = / <* «WI = 7 *L ^W . (B .64) 



u J n p' — p — te 

4 m 2 —oo 

Therefore, if considered as a function of p with s = 4 (m 2 + p 2 ), the chiral 
moment V^J R (p) is analytic in the upper complex half plane. We note that 
the full vertex function V v ^{s) on the other hand is not analytic in the upper 
complex half plane if considered as a function of p. It shows branch points at 
p = ±im which are systematically eliminated performing the 1/m expansion. 



Start with the first term in (B.62). The result can immediately be deduced 
from the simple identity: 



3 In 2 (l-2i— ) = - arctan f 2 — ) In ( 1 + 4^-) . (B.65) 

The second term follows with the somewhat more involved identity: 



^ . , mm n rnt ^ ( /l ip\ (l m^ \ it 2 
XR mr b) = - r ^-z±9t[lA 2 U + — -2Li 2 *— + 

47r 8p d \ \2 mnJ \2m n — ipJ 12 

2 / m n — 2ip\ 1 2 / 1 % p 



_ ln2 -^ -^ + _ ln ^ _ _ J21 . (B.66) 
\2m n -2ip) 2 \2 m^J J v ; 

For both terms we succeed in writing the spectral function as the imaginary 
part of an analytic function with branch cuts only on the lower complex half 
plane. Thus it is straightforward to reconstruct the full analytic function. We 
collect all terms and display our result 9 : 



T/ (2) / x m 2 ml ml ( (I ip\ (\ m w \ . (I 

- In 2 f™*~ 2 *M - In f i - ^) In 2^ . (B.67) 

\2m n -2ip \2 m w J I v ; 



B.8 Three-loop bubble 

We turn to the three- loop bubble 

J 7rw ( s ) = -ltr|^ i r,r 275 C- 1 ^(P + /)V^(P + /,^-/)^(^-/). 

(B.68) 

which is defined in terms of the full off-shell two-loop vertex of (B.60). We 
introduce the loop-subtracted bubble function, J nn (s), with 



J wn (s) = J^(s) - 2 c(s) J(s) J K (s) + c 2 {s) J 3 {s) + ■■■ (B.69) 

where we dropped terms of suppressed chiral orders induced by one and two- 
pion production. The loop J niT (s) is decomposed into its chiral moments: 



Note the identity: Li 2 (l) = f| - \ In 2 2 
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^7T7r\ $ J 



E 

n=3 



Jil\p) 



(B.70) 



where J^(p) carries chiral power n. The spectral function is expressed in 
terms of the renormalized one and two-loop vertex functions V n ^(s) and 

^rw,fl(«) 



3 J_(s) = 2 (9 J(s)) (ft K^(s)) + (9 J(s)) (K iR ( S ) \/; jR (s)) . (B.71) 
For the leading moment Jl 3 J(p) we find in analogy to (B.35) the relation: 



Ji$(p)-J!8(o)=i^(2vV(p) + 



2tt 



V. 



(1) 



>) 



(B.72) 



r(l)|2/ 



Kl)^M2 



where we introduced the proper analytic continuation |lv#|e(p) of |V^(p)p 
with: 



V$(p)K ( X'*(p) = B 



V, 



(1) 



7T,R 



!(P) 



= ]4?v r U J - Ll2 U + ^ 

+ In (2) In ( 1 - 2 — ) + ^ In (16) ) . 
V 77W m*- ' / 



(B.73) 



We collect all terms and present the leading chiral order contribution 



" W ™ V ; (4tt) 3 p \ \2 



1 m, 



7TT.7T — ip 



ip 



m-n 



In (16) + In 



2 / m n — lip 



T\- 



2m v — lip J 6 
Note that the loop function J^J(p) is convergent with 



(B.74) 



JS > ( o)./'^Mf<£U n(16 ) 



71 p^ 



m ml 



(4vr) 



(B.75) 



The integral (B.75) follows from (B.74) and the observation J^(oo) = 0. 
The leading order renormalized three-loop function is J nir r(p, z) = J^J(p) — 
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C NNLO terms in CE scheme 



In this appendix we systematically collect all terms of chiral order one in 
CE scheme. We emphasize that again explicit calculation shows that no in- 
dependence survives at this order. The auxiliary coupling function gj^ receives 
the contributions 



9h{s, m w ; c) = g^ s (s, m w ; c) - J(-i z ) - g\ A 2 J w>s (-i z ) 

-^A 4 J_, s (-^o)-(c-i)%iA 6 J 4 (-^ ) + --- (C.l) 

where we included the anomalous contribution from J n7Tn ^s according to (73). 
It is convenient to resolve the structure of the subtraction constants J<x,s{~i z o) 
and J-k-k s{— izo) with: 



9a A 4 J *,s(-i zo) = (c - I) ^j{zq) + g\ A 2 (kj^Zq) + ^l{z )) , 
9a^ 6 J™,s(-i zo) = (c - |) Kj(z ) ((c - |) k] (z ) + 2 g\ A 2 Rj^(z )) 

+ g A \ 2 Kj^(z ) (C.2) 

where we introduce Rj(zo) = g A \ 2 J(—iz ), Kj^{zq) = g A \ 2 J n (—izo) and 
Kj^n(zo) = g\ A 4 J 7r7 r(— izo). For later convenience it is appropriate to gener- 
alize here our notation and introduce further dimensionless coefficients R Jn k (z) 
and K^ nk {z): 



^Az) = g A k \ 2k (j£)' ±d l p J^ R {-iz,z) , 

K^(z) = g A k \ 2k Q£) ^d l p V„ k>R (-iz) . (C.3) 

The auxiliary coupling gn does show an explicit c-dependence at chiral or- 
der Q 3 . However such terms are canceled by a contribution to the scattering 
amplitude with an explicit c-dependence 



The cancellation follows from the intermediate result 



g A \ 2 g- H \s)-2Rj{zo) (c-\) g- H \s) = g 2 A \ 2 g~ R \p) + O (Q') (C.5) 



which considers the relevant combination. Again our result (C.5) relies on 
the leading order expression (92). We decompose the amplitude in its chiral 
moments T[i So ](p) = X^L-i T^Jp) and display our complete collection of 
terms with chiral power one: 



1 0J Jr(p,z) 27r Jr(p,z) 2ttJ^(p,2;) 

Qr 2 (p) /- 9\ 2 fn 



JUp,z) ^ «' ^ 2 J 3 R ( P ,z) 

A? rr? — iR^PL _ „4 x 4 ( c, *W0) J**,r(P, Zq) 

+ H2m -2n J R (p, z) 19aX \ l J R (p, z ) ' J%(p, z ) 

_ 4 2 A 2 fl^O) ( VnAP) _ J*,r(P,Zq) \ 

A Jr(p,z) \Jr(p,z ) Jr{p,z ) J 

Q 4 A 4 \2 

- Z-Fzf ^{ J r(P' z o) K,r(p) - J tt,r(p, zo) 

J R\Pi Z 0) 
,4 \4 



+ g\X i [B 7T7T (u,t) + B 7T7T (t,u)) . (C.6) 

Our result (C.6) involves the previously given renormalized one-loop ver- 
tex V^Ap) ~ Q an d two loop bubble J^APi z o) ~ Q 2 - It further requires 
the renormalized two-loop vertex V^^Ap) ~ Q 2 an d the three-loop bubble 
Jtttt,r(p, zo) ~ Q 3 - Again all higher order pole terms must cancel. This re- 
quirement confirms C2 = C2 an d determines the hitherto unknown coefficients 



n _ KjM z o) 1 7 K %^) rn - m 

m% 777.Tr 2 7T 2 7T 

in terms of the known parameters (2 and Kj^(zo), Kj j7r7r (zo) 10 . 

We make the analytic structure of our result (C.6) more explicit by expressing 
our result in terms of residual functions Rj^h (p, z) and Ry ' k (p, z) which are 
defined 



Tl-l 

g 2 / * 2k J^R(p, z)=g\ k \ 2k J n R (P, z) R<$.(p, z) + Y, *%*{z) J l R (p, 

1=0 

9a * 2k V nk Ap)=9 A k A 2 ' JR(p, z) R { ;l(p, z) 



10 Note that here we suppress a finite renormalization of Q ~ * Co \X ~ ^ z o C2 , 
This renormalization leads effectively to g^ ~ Q (p 2 + 2;2 )- 



ra-1 

+E4/W4(^) (c.8) 

in terms of the coefiicients ftj„.fc(z) and R Vnk (z) of (C.3). This leads to the 
central result of this appendix 



r, ( .>) = - -^ + - Cf ' (2 4°iM - 2 *£(») - 3 Cf ' « 

i jo) m p + « z m v 

- — (d 0) ) 2 ^ + 2 ^A 4 ^(p) 

,2 ™/d(3)/„ >>_ p(2) 



+ 4 <?i A 2 tf ' < f- ^>, *>) - R { vUPi *. 
- 4^ A 2 (i?g(p, zo) - <lb, «,)) d 0) (zo +ip] 



+ 2 <£ A 4 (* -ip)^ «>> *>) - 4 2 i(P, ^o)) 2 • (C.9) 

where we decomposed the residuum w = J2T=o w m ^° its chiral moments. In 
(C.9) we encounter the reduced functions <L R%„ and R™, R% defined 
according to (100). Note that all terms in (C.9) are regular at p = — iz and 
p = —izQ except the first one. The contribution from the pion nucleon box 
diagram sits in R nix {p) (see (80)). With (C.9) we arrive at the desired repre- 
sentation of the scattering amplitude. The chiral correction terms of order one 
renormalize the pole residuum w and add further smooth remainder terms. 
The pole residuum is: 



l-w = R%(z )-2 4l(z ) + 2^ 0) z 

+r%M) - 2«i?L(*b) - (4°i(*o) - 4>o)) : 
+ 4(4 t(^o)-4:i(-o)-d 0) ^o)d 0) ^o 

~V)U\ rJ 2 )f,\\ m X 2 )™2 , n (rfl 



2 (R^(z ) - «$(«,)} ^ C K + O (Q*) . (CIO) 



We close this appendix with a more detailed analysis of the low energy be- 
havior of (C.9). A systematic zo/m^ expansion is worked out. In particular we 
show results for the low energy coefiicients b e g and b^ Q ■: 



, (o) _ _ >(o) _ rn 2 2 / (i) (2 ) 



^ = -C^ - f^ 9 a A 2 (1 - 4 R^(0, z ) + 2 i2£(0, ^ 
-d 0) (2 4°lU)-2^i(,o)-3d 0) .o 
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-9U 2 tf ] rnl^(Rfl(0,z )-RW(0,Zo: 

2 71 K 



,2 



A 2 (4i(-o) - «£(«,) - d° } *b)^ (4>, zo) - <(0, ,o 



9A v , 

- <?! a 4 ^ (4 2 L(o,^o) - 2 <L(o,z ; 

777 / 

$ = ~9\ A 2 £ H $,) Ki(0, *b) - 2 <t(0, z )) + (d 0) ) 2 



-«a A 2 ^ H $,) (i2JS(0, z ) - 2 it^(0, z )J + (C: 

- 9 2 a A 2 fj Ci 2) ™ 2 (-< P ) (i?g(0, zo) - 4 2 t(0, ^o 
+ ^A 2 -d 0) (4 2 i(0^o)-<(0,, 



Ti 

+ g 2 A X 2 - Cf z (-» d p ) (i2g> (0, z ) " 4>, ^ 

-^A 2 (4°t(^ ) - «gl(*b))^ (-^ p ) (4 2 i(o,^o) -4>,z ; 



<?1 A 4 — (-< $,) (7^(0, z ) - 2 7^(0, zo) 
<?1 A 4 ^ (-< P ) «(0, zo) - 4>, *bf 2 



2 , .. , 2 



47T 

+ ^A 4 ^« ) 7r (0,z )-4 2 i(0^o 

- si A 4 ^ H «9 P ) 7^(0) + O (Q 2 ) (C.ll) 



introduced when discussing the structure of the scattering amplitude in (27). 
First the coefficients K Jir k (z) and Ky k ( z ) are systematically expanded in the 
small ratio z/m n with the leading moments 



foi / n 9 ^ 9 m m ?r / -2 4 z 2 2 z 3 „„ / z 4 N 

4°t(z) = ^ 2 A 2 — -^ 1 + — + - — + -— + — 



4 7r y m^r 3 m^ m 3 \ mi 

fil/ n 2 , 2 /I 4z 3z 2 32z 3 m f z 4 ' 

« z = g\ A 2 - + + + +o( — 

\ 2 3 m T mi 5 m; V m* 



^ U) "^ A (4tt) 2 [n [ml V + U {A'm„ 



4^)=^A 2 ^ 2 + ^ + ^ + ^f + O ^ 
' 4tt V m^ ml ml ml, 



(2) , x ,,,/ 4 z 12 z 2 32 z 3 m ( z 4 ' 



m^ 7774 m% \ mi 
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,(0) / >»_„4 X 4 m2m2 ~ 



«^W=^A 4 — -f In (2) + 

[Airy \ 2 m n 



RSi.W = 9 jA^Jl + 2 1n(2) + (f + f ln(2))e- 




2 2 




+ (l4+121n(2)) S j + 0(^)). (C.12) 

Next we provide the leading moments of the residual functions 



4l(p,z 






R%{p,z 




n(2) 

n .J,TT7T 



.(C.13) 
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Finally note the low energy characteristics of the s-wave projected box function 



m m n / 1 % p ( p 2 

Att I 4 2m T [m 2 



ft.(p) = ^l7 + ^- + olS- ■ < C14 ' 



Applying the our results (C.12) and (C.13) we derive 



l_„, = ^ si A'(2*. + ^4 + o(4)) +2 C f* 

4vt \ mj,- 3 7774 \m%J J 

+ (^)'((-T + i h »); 
+ (-^ +8hW )i +0 (i)) 



ram^ ^ j 3zq | g /j£_^ ^ A ( ) 
4 vt \ m^ \ m 2 . 



-MJV1 + - + OR (i% 



4vt y dm T m; \m*J J 



• 



:-! 



4vt \ 3m ff m; \m*J J 

+ ("'^>)lHti)) 

2 A 2 m? "^ (2) /8 | 16 z | 288^ 2 | g /4\\ 
4 7r 1^3 rn^ 5 m 2 , \ m w// 

+ ^A^cf(2 + ^ + 5ii + (4)) 

4 7r \ m n Jm; \m*J J 



J 

/ 1 fi -V. / -V 



>2 = -^Ua 2 (i + ^ + o(4 

7r ml \ 5 77V V 7771 



4,4 "i 2 / /547 64 , . A z 

+ 9a A t r^ 1 + In (2) — 

yA 4tt 2 I V 30 5 V V m^ 



/2741 64 

+ 



V 45 



>v)i +o (S)) 



4 7r y om^ 7774 ym^. 







{d 0) ) 2 +o( Q >). 





(C.15) 
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